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Preface 


The  purpose  of  this  study  was  to  determine  if  there  existed  an 
analytical  solution  for  the  three-body  atmospheric  entry  equations.  I 
decided  to  continue  the  excellent  work  by  Drs.  Vinh,  Busemann,  and  Culp 
and  use  the  Method  of  Matched  Asymptotic  Expansions  to  determine  zero 
and  first  order  solutions  to  the  above  equations.  By  using  the  computerized 
symbolic  manipulator  Mathematica™,  the  algebraic  manipulations  are 
reduced  to  a  manageable  level,  so  a  more  thorough  analysis  can  be 
performed.  By  examining  trends  in  these  higher  ordered  solutions,  I  hoped 
a  complete  analytic  solution,  not  evident  from  traditional  analytic  means, 
would  become  evident. 

I  would  like  to  thank  many  people  whose  work  aided  my  study. 

First,  I  would  like  to  recognize  my  faculty  advisor,  Capt  Rodney  Bain. 
Without  his  keen  mathematical  insight ,  continuing  motivation,  and 
unfailing  confidence,  this  study  would  not  be  possible.  I  am  greatly 
indebted  to  the  work  of  my  predecessor,  Harry  ICarasopoulos.  Much  of  the 
groundwork  of  this  study  is  based  on  Harry  's  master's  thesis  at  AFIT.  I 
would  also  like  to  thank  Drs.  Vinh,  Busemann,  and  Culp,  whose  work  in 
analytical  flight  mechanics  greatly  aided  me  in  my  study,  and  also  Dr. 
Stephen  Wolfram,  whose  Mathematica™  code  saved  me  many  late  evenings 
from  algebraic  tedium.  Finally  I  would  like  to  thank  my  parents  for  their 
invariable  motivation  and  support  throughout  my  education.  Without  them, 
this  study  would  not  be  possible. 


Ted  Masternak 
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y  Flight  Path  Angle  (deg  or  rad) 
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1  Inner  Expansion 

c  Composite  Expansion 

0  Outer  Expansion 
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Inertial  Time  Derivative 
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Abstract 

Although  numerical  techniques  exist  to  obtain  solutions  to  highly  non¬ 
linear  and  highly  coupled  systems,  the  trends  and  subtleties  of  the  solution 
are  frequently  lost  in  the  volume  and  form  of  tabular  and  graphical  data  in 
covering  a  wide  range  of  initial  conditions.  By  deriving  an  approximate, 
analytical  solution,  relationships  between  dependent  parameters  are 
discernable.  Also,  the  derived  solution  is  easily  applied  to  any  new  set  of 
initial  conditions  or  can  be  modified  to  incorporate  slightly  different 
equations  of  motion.  This  study  presents  an  analytical  investigation  of  the 
three-dimensional  equations  of  motion  for  lifting  entry  into  a  planetary 
atmosphere. 

In  this  study,  the  equations  of  motion  for  lifting  entry  into  a  planetary 
atmosphere  are  derived.  A  non-rotating,  spherical  planet  is  assumed,  as  is  a 
non -rotating,  strictly  exponential  atmospheric  model.  The  derived 
equations  of  motion  are  transformed  to  a  variable  set  relating  the  classical 
orbital  elements  to  the  vehicle's  altitude.  Solutions  to  the  resulting  five  non¬ 
linear,  coupled,  first  order,  ordinary  differential  equations  are  obtained  by 
using  the  Method  of  Matched  Asymptotic  Expansions  and  a  computerized 
symbolic  manipulator,  which  performs  the  detailed  algebraic  computations. 
By  using  the  planetary  scale  height-mean  equatorial  radius  (PSHMER) 
product  as  a  small  parameter,  both  zero  and  first  order  expansions  to  the 
equations  of  motion  are  obtained. 

viii 


It  is  demonstrated  the  analytical  solution  agrees  with  results  obtained 
from  numerical  integration  of  the  equations  of  motion.  Due  to 
approiimations  made  in  the  solutions  of  the  first  order  inner  expansions,  the 
analytical  solution  slightly  deviates  from  the  numerical  solution  at  tow 
vehicle  altitudes.  The  two  solutions  are  compared  further  and  the  validity 
of  the  analytical  solution  is  examined. 
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USING  THE  METHOD  OF  MATCHED  ASYMPTOTIC  EXPANSIONS. 
ANALYTICALLY  INVESTIGATE  THE  THREE-DIMENSIONAL. 
ATMOSPHERIC  ENTRY  PROBLEM 


I.  Introduction 

Since  the  recent  explosion  of  readily  available  computing  power, 
numerical  solutions  of  complex  systems  described  by  non-linear  phenomena 
have  become  ubiquitous.  In  using  these  numerical  techniques,  the  trends 
and  subtleties  of  the  original  system  often  become  lost  in  the  results  of 
numerical  analysis.  By  deriving  a  simple  and  relatively  accurate  analytical 
solution  to  a  complex  physical  system,  a  readily  available  analysis  becomes 
available  which  retains  the  subtleties  of  the  original  system 

The  equations  of  motion  for  lifting  entry  into  a  planetary  atmosphere 
are  highly  coupled  and  non-linear.  Instead  of  using  a  numerical  technique 
to  solve  the  equations  of  motion,  an  approximate  analytical  method  can  be 
used,  where  the  trends  and  patterns  of  the  original  system  are  preserved 
and  expressed  in  a  solution  which  is  simple,  accurate  and  practical. 
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In  this  effort,  the  three-dimensional  exact  equations  of  motion  for 
lifting  entry  into  a  non-rotating  planet  are  developed  to  first  order 
accuracy.  Solutions  to  the  equations  of  motion  are  developed  by 
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considering  the  atmosphere  as  a  boundary  layer  which  perturbs  the 
vehicle  s  motion  as  is  approaches  the  planet's  surface.  The  Method  of 
Matched  Asymptotic  Expansions  is  used  to  develop  zero  and  first  order 
solutions  valid  throughout  the  flight  trajectory  of  the  entry  body.  Due  to 
coordinate  singularities,  this  study  is  limited  to  non-potar  and  non- 
equatoria!  entry  trajectories. 

Assumptions 

In  this  study,  the  non-rotating  planet  is  modeled  as  a  spherical  body 
possessing  a  inverse  square  law  gravitational  field.  A  strictly  exponential, 
non-rotating  atmosphere  is  assumed  with  its  density  as  a  function  of  the 
radial  distance  form  the  planet  s  surface.  This  investigation  assumes  the 
vehicle  s  trajectory  is  influenced  only  by  aerodynamic  forces  and  the 
planet  s  gravity.  Other  perturbing  forces,  such  as  solar/lunar  gravitational 
forces  and  other  celestial  perturbations,  are  considered  negligible.  The 
vehicle  s  lift-to-drag  ratio  and  ballistic  coefficient  are  assumed  constant  and 
prescribed. 

Methodology 

In  Section  II,  the  equations  of  motion  for  three-dimensional,  lifting 
entry  for  a  spherical,  non-rotating  planet  are  derived  from  basic  kinematic 
and  force  relationships.  In  Section  III,  these  equations  of  motion  undergo 
coordinate  transformations  to  express  the  independent  and  dependent 
variables  as  convenient,  non-dimensional  orbital  parameters.  The  first 
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order  accurate  analytical  solutions  are  derived  in  Section  IV  using  the 
Method  of  Matched  Asymptotic  Expansions.  The  five  coupled,  non-linear 
ordinary  differential  equations  of  motion  are  expanded  to  two  separate,  but 
overlapping  domains:  the  outer,  Keplerian,  domain  and  the  inner, 
atmospheric,  domain.  Each  of  these  expansions  is  solved  independently  to 
create  two  sets  of  solutions  to  the  original  equations  of  motion.  Finally  the 
two  solution  sets  are  blended  together  to  generate  one  solution  valid  over 
the  entire  vehicle  trajectory.  The  validity  of  these  analytical  solutions  is 
examined  in  Section  V.  It  is  shown  the  analytical  solutions  encounter 
singularities  near  polar  orbits  and  slightly  underestimate  aerodynamic 
turning  at  low  altitudes.  Section  VI  summarizes  the  study's  findings  and 
presents  recommendations  for  further  study. 
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II.  Derivation  of  the  Equations  of  Motion 

Introduction 

In  this  section,  the  original  equations  of  motion  are  derived  for  three- 
dimensional  lifting  entry  into  a  non-rotating  planetary  atmosphere.  The 
non-thrusting,  lifting  vehicle  is  modeled  as  a  point  mass  in  a  three- 
dimensional  coordinate  space.  The  vehicle's  orbit  is  assumed  to  be  initially 
described  by  Keplerian  or  two-body  motion.  Aerodynamic  forces  are 
assumed  to  be  the  only  perturbations  acting  on  the  entry  vehicle. 

For  the  sake  of  brevity,  the  derivations  in  this  section  are  the 
abbreviated  versions  of  much  more  detailed  work  presented  in  other  studies 
(Vinh  and  others,  1980:20-28)  and  (Karasopoulos,  1988:9-36).  Thorough 
descriptions  of  the  derivations  are  presented  in  the  two  excellent  references 
listed  above. 

Coordinate  Systems 

Figure  1  defines  the  planetary  coordinate  systems  used  in  this  study. 
In  addition  to  coordinate  systems  referencing  the  vehicle  to  the  planet, 
there  is  also  present  a  coordinate  system  centered  on  the  vehicle  and 
relating  its  orientation  to  the  planet  below.  Listed  below  are  the  coordinate 
system,  their  respective  coordinate  and  unit  vector  representations  and  a 
brief  description  of  each  (Vinh  and  others,  1980:22-24). 
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Figure  1.  Planet  Centered,  Rotating  Coordinate  System 


Reference  Unit  Rotation 

Coordinates  Vectors  Rates  Description 

A  «  * 

OX,YjZ|  I J  1C  None  Inertial  frame.  Inertial  reference  frame 

whose  center  is  coincident  with  the 
center  of  the  spherical  planet  s 
gravitational  field.  The  OX,Y,  plane  is 
in  the  equatorial  plane  and  the  Z,  axis 
completes  the  right-handed  system. 

OXYZ  ipjpkp  oo  kp  Planet  frame.  Non -inertial  reference 

frame  whose  center  is  coincident  with 
the  center  of  the  spherical  planet  s 
gravitational  Held  and  fixed  with 
respect  to  the  planet.  The  OXY  plane  is 
in  the  equatorial  plane  and  the  Z  axis 
completes  the  right-handed  system. 
Thus,  the  coordinate  system  rotates  at 
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the  same  rate  as  the  planet,  w,  about  the 
kP  axis. 


*  »  *  *  * 

Oxyz  iB  jB  kB  6kp,  4>jB  Body  Frame.  Non-inertial  reference 

frame  whose  center  is  located  at  the 
vehicle.  The  x-axis  is  along  the  position 
vector  from  point  0  to  the  vehicle,  the 
y-axis  is  in  the  equatorial  plane  and 
orthogonal  to  the  x-axis  and  the  z-axis 
completes  the  right-handed  system. 


Ox  y'z'  >wJw^w  ^iwYkfl  Wind  frame.  Non-inertial  reference 

frame  whose  center  is  located  at  the 
vehicle.  The  x-axis  is  along  the  lift 
vector  (o  -  0)  from  point  0  to  the 
vehicle,  the  y-axis  is  along  the  drag 
vector  (a  -  0)  and  orthogonal  to  the  x- 
axis  and  the  z-axis  completes  the  right- 
handed  system. 


Applying  the  coordinate  transformations  to  the  above  systems  results 
in  two  coordinate  transformation  matrices.  Combined,  the  two 
transformations  relate  the  vehicle's  coordinate  system  in  terms  of  the 
planet's  rotating  coordinate  system. 


X 

'  cos(0)  cos($)  -sin(8)  -sin(4>)  cos(0) ' 

x' 

y 

- 

sin(0)  cos($)  cos(0)  -sin(0)  sin(4>) 

y 

(2.1) 

z 

p 

sin($)  0  cos(4>) 

z. 

B 

i 

cos(y)  sin(y)  0 

x' 

y 

- 

-sin(y)  cos(<j> )  cos(y)  cos(<J>)  -sin(<|>) 

y 

(2.2) 

z 

B 

.  -sin(y)  sin((j>)  cos(y)  sin(<|>)  cos(i|>) . 

z. 

w 

Kinematic  Equations 
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The  vehicle’s  position  vector  is  given  as  r  -  riB.  Differentiating  this 
position  vector  to  find  the  inertial  time  derivative  of  the  position  vector  uses 

^d?  Bn?  Di 

the  vector  differentiation  relationship  ■jjj-  -  —  +  co  x  where 

-,B1  •  . 

u>  -  ((co  ♦  0)  sin(4>))  iB  -  4>  jB  ♦  ((u>  ♦  0)  cos(4>))  kB. 


d l 

dt 


(dr) :  (  /xi  (d0 

dt  rcos(4>)  j+w 


\  * 


Jb 


me 

dt  kB 


(2.3) 


The  vehicle  s  local  velocity  with  respect  to  the  Oxyz  reference  frame 
uses  Eq  (2.2)  and  is  given  as 

^  *  vJw  "  (Vsin(y))iB  ♦  (Vcos(y)cos(<J>))jB  ♦  (Vcos(y)sin(<j>))kB  (2.4) 
Equating  Eqs  (2.3)  and  (2.4)  gives  the  kinematic  equations  as 

dt- Vsin<y)  (2.5) 


d0  V  cos(v)  cos(d) ) 

dt  ’  r  cos(4>)  " 00  (2  6) 


d<t>  V  cos(y)  sin(d)) 

dt "  r -  (2.7) 


Force  Equations 

Since  the  acceleration  and  the  force  on  a  vehicle  are  related,  the 
acceleration  of  the  vehicle  is  calculated  in  deriving  the  force  equations. 
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Taking  the  inertial  derivative  of  Eq  (2.4),  substituting  in  Eqs  (2.6)  and  (2.7) 
and  solving  for  the  derivatives  dV/dt,  dy/dt  and  d<i>/dt  gives 


dV  Ft 

dt-  m  *sin(y> 


*oj2r  cos(<j>)  (sin(y)  cos(A)  -  cos(y)  sin(<t>)  sin(i|i))  (2.8) 


dy  Fn  cos(o)  V2 

V  dt  *  m - -  g  cos(y)  +  ~  cos(y)  +  2ooVcos(4>)  cos(<J> ) 

+  u>2r  cos(4»)  (cos(y)  cos(<t>)  +  sin(y)  sin(4»)  sin(<|>))  (2.9) 


d(b  Fn  sin(a)  v2  co2r 

v  dt  “  m  cos(y) "  T  cos(y)  cos^  >  tan(<i>) '  cosjy)  sin(<}>)  cos(<t>)  sin(^ ] 

+  2coV(tan(y)  cos(<J>)  sin(i|> )  -  sin(<J>))  (2.10) 

Fn  is  defined  as  the  lift  force,  FT  is  defined  as  the  aerodynamic  and 

propulsive  forces  along  the  velocity  vector  and  o  is  called  the  bank  angle, 
which  rotates  the  FN  vector  out  of  the  local  vertical  plane. 


Assumptions 

The  above  three  equations  are  known  as  the  force  equations.  Since 
this  study  assumes  a  non-thrusting,  lifting  vehicle,  the  vehicle  s  thrust  is 
zero  or  T  -  0  and  by  the  definitions  of  FT  and  FN  ,  FT  -  -D  and  FN  -  L. 

Since  this  study  assumes  the  planet  and  its  atmosphere  are  not 
rotating,  co-  0.  This  assumption  is  commonly  used  in  entry  flight 
mechanics,  where  the  analysis  is  primarily  concerned  with  the  variations  in 


the  velocity  and  altitude  of  the  entry  vehicle  in  the  portion  of  the  trajectory 
where  high  deceleration  develops  (Vinh  and  others,  1980:27).  In  his  study. 
Karasopoulos  showed  exclusion  of  a  rotating  planet/atmosphere  results  in 
an  analysis  not  valid  for  some  entry  trajectories,  especially  where  the 
vehicle  undergoes  a  long  atmospheric  entry  period,  as  in  very  shallow  entry 
trajectories.  This  study  assumes  the  non-rotating  equations  of  motion  are 
valid  and  applies  them  to  trajectories  where  they  are  accurate. 

Substituting  these  assumptions  in  the  above  equations  result  in  the 
modified  kinematic  and  force  equations. 

^-Vsin(y)  (2.11) 


d6  V  cos(y)  cos(4>) 
dt  ’  r  cos(4») 


d£  V  cos(y)  sin(^) 
dt  "  r 


V 


V 


dV  D  ,  > 

dT--m-gs,n,y) 


dy  L  cos(g)  ,  v  V* 
dt  "  m 


-  g  cos(y)  ♦  “jrcos(y) 


d<l>  L  sin(o)  V2 

4\  "  m"cos(y) '  T  cos(y)  cosi(b ]  tan(<t>) 


(2.14) 


(2.15) 


(2.16) 


Summary 
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This  section  briefly  derives  the  equations  of  motion  (the  equations  for 
flight  over  a  non -rotating  spherical  planet)  for  a  vehicle  considered  as  a 
point  mass  traveling  within  a  planetary  atmosphere.  The  derived  equations 
are  repeated  below  for  completeness. 


57  -  v  sin(y) 

(2.11) 

d0  V  cos(y)  cos(4>) 
dt  ’  r  cos($) 

(2.12) 

d<t>  V  cos(y)  sin(<l>) 

dt  "  r 

(2.13) 

dV  D  ... 

^-----gs.n(y) 

(2.14) 

dy  Lcoslo)  .  .  V2 

dt  *  m  -  gcos(y)  +  -jrcos(y) 

(2.15) 

d<l>  L  sin(o)  V2 . 

dt  •  m  cosly) '  r  cos(*)  cos(*>  “"<♦> 

(2.16) 
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III. 


The  equations  of  motion  for  three-dimensional,  lifting  entry  into  a 
non-rotating  planetary  atmosphere  were  derived  in  Section  II.  In  this 
section,  assumptions  and  approximations  are  defined.  Coordinate 
transformations  are  undertaken  to  convert  the  equations  of  motion  into  a  set 
of  coupled  ordinary  differential  equations  with  convenient  dependent  and 
independent  variables  The  equations  of  motion  derived  in  Section  II  are 


given  as 


■If  -  V  sin(y) 


(3.1) 


dO  V  cos(y)cos(<l>) 
dt  "  r  cos($) 


(3.2) 


d£  V  cos(y)  sin(d>) 
dt  r 


(3.3) 


dy  D  ... 


(3.4) 


dy  L  cos(q) 
dt  "  m 


g  cos(y)  ♦  —  cos(y) 


(3.5) 


(3.6) 
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Assumptions 


Spherical  Planet/Atmosohere.  Although  planets  are  usually  oblate 
ellipsoids,  the  spherical  planet  assumption  is  common  in  entry  flight 
mechanics  analysis,  since  the  ellipticity  of  the  planets  are  of  the  order  I  O'2 
to  10'4  (Vinh  and  others,  1980:3).  Thus,  the  error  induced  by  this 
approximation  is  insignificant. 

Following  the  spherical  planet  assumption,  the  planetary  atmosphere 
is  assumed  to  be  a  sphere.  In  reality,  an  atmosphere  is  approximately  an 
oblate  ellipsoid  with  minor  deviations  due  to  solar  storms  and  uneven 
heating  of  the  planetary  atmosphere.  Generally,  these  effects  are  present  at 
altitudes  above  250  kilometers,  where  aerodynamic  forces  are  insignificant 
except  in  the  slow  orbital  decay  of  satellites  (Vinh  and  others,  1980.2). 


Gravitational  Model.  Since  the  planet  is  being  modeled  as  a  sphere, 
its  gravitational  Held  follows  Newton  s  inverse  square  law  and  is  given  as 


g(h)  •  g. 


(3.7) 


Atmospheric  Density  Model.  The  planet  s  atmosphere  is  assumed  to 
be  a  non-rotating  sphere  fixed  with  respect  to  the  planet.  From  the  state 
equation  for  a  gas  and  the  hydrostatic  equation,  the  planetary  atmosphere  is 
modeled  by  the  equation 


d£ 

P 


-fj  dr 


(3.8) 
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where  $  is  the  local  atmospheric  density  and  r  is  the  radial  distance  from 
the  center  of  the  planet.  The  1/fJ  term  is  called  the  scale  height  and 
describes  the  size  of  an  altitude  region  where  the  atmosphere  is  accurately 
modeled  by  an  exponential  relation. 

As  in  many  other  analytical  entry  flight  mechanics  studies,  it  is 
assumed  the  quantity  £r  is  a  constant  for  a  planetary  atmosphere  (Vinh  and 
others,  1980:5).  For  most  planets,  this  value  is  usually  of  order  1000.  For 
earth,  it  has  an  average  value  of  about  900.  Thus,  solving  Eq  (3.8)  with  this 
assumption  yields 


P-P5 


(3.9) 


This  study  assumes  the  planetary  atmosphere  is  accurately  modeled 
by  a  exponential  atmospheric  model  given  as 


-fry 

p-p,e 


(3.10) 


where  y  is  the  altitude  above  the  planet  surface.  By  definition,  h  -  y/r, 

Thus  (ty  -  (3  hr,  -  h/c  or 

-h/e 

P-P,  e  (3.11) 


where 


The  value  of  €  is  on  the  order  of  1/1000  and  thus  is  a  small  number.  Using 
this  variable  as  a  small  parameter  allows  the  use  of  Matched  Asymptotic 
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Expansions  to  be  applied  to  derive  an  analytical  solution  to  the  equations  of 
motion. 

Aerodynamic  Forces.  The  lift  and  drag  parameters  for  the  entry 
vehicle  are  given  by  the  following  relationships,  which  incorporate  the 
exponential  atmospheric  density  relation  determined  above. 


pse'h/<  V2S 


(3.13) 


evfs  p,e~l,/e  v2s 

D  -  CD  2  '-'D  2 


(3.14) 


Although  CL  and  C0  are  functions  of  the  angle  of  attack,  Mach  number 
and  other  flight  parameters,  this  study  assumes  they  are  constant  and 
prescribed  at  the  beginning  of  the  entry  trajectory.  This  is  a  common 
assumption  when  dealing  with  hypersonic  flight  mechanics  (Vinh  and 
others,  1980:101). 

Substitution  of  these  assumptions  into  the  equations  of  motion  given 
by  Eqs(3.1)-(3.6)  gives 

^-Vsin(y)  (3.15) 


d0  V  cos(v)  cosM) 
dt "  r  cos($) 


d<j>  V  cos(y)  sin(d>) 

dt "  r 


(3.17) 


dV  c„  p,e-h/‘  V2S  {r%\2 

df  25 - "  *•  (r  J 


sin(y) 


(3.18) 
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(3.19) 


h/<V2Scos(a)  (r.\Z  y2 

~2^ - "g*  r  cos(y) +  ycos(y) 


CL  Pse'h/c  V2S  sin(c) 


dt 


2m  C0S(y) - -  —  cos(y)  cos(4>)  tan(4>)  (3.20) 


Transformation  of  the  Independent  Variable 


As  given  in  Eqs  (3.1M3.5),  the  equations  of  motion  are  in  terms  of  the 
dimensional  variable  time,  t.  For  atmospheric  entry,  it  is  more  convenient 
to  relate  the  motion  of  the  vehicle  in  terms  of  the  orbital  altitude,  r,  since 
the  two  major  forces  controlling  the  entry  trajectory,  gravity  and  lift/drag, 
are  functions  of  r.  To  change  the  integration  variable  from  t  to  r.  the  chain 
rule  for  differentiation  is  used.  Thus,  using  dr/dt  relationship  from  Eq  (3.1 ), 
differentiation  with  respect  to  r  is  defined  as 


A.  d  dt  _  1  d 

dr  "  dt  dr  "  V  sin(y)  dt 


(3.21) 


This  transformation  reduces  the  number  of  equations  of  motion  from  six  to 
five  since  Eq  (3.1 )  is  incorporated  into  the  other  five  equations  by  the  above 
use  of  the  chain  rule. 

In  atmospheric  entry  analysis,  a  dimensionless  altitude,  h,  is 
commonly  used  as  the  independent  variable  of  integration  for  the  equations 
of  motion,  h'  is  defined  as  (Vinh  and  others,  1980:256) 

h  -  f  (3.22) 
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Since  r  is  currently  the  integration  variable  for  the  equations  of 
motion,  a  transformation  relationship  is  needed  to  convert  r  to  h.  By 
definition 


r  -  rs  ♦  y  -  rs  ( 1  +  h) 


(3.23) 


Differentiation  the  above  equation  gives 


dh”r* 


Again  using  the  chain  rule  for  differentiation  yields 

d  d  dr 
dh  ’  dr  dh  “  r*  dr 


(3.24) 


Applying  the  above  definition  for  r  as  well  as  the  differentiation 
transformation  for  d/dr  yields 

d  d  r.  d 


dh  5  dr  ’  V  sin(y)  dt 


(3.25) 


Applying  this  transformation,  and  the  definition  of  r,  to  Eqs  (3.15)- 
(3.20)  gives  the  modified  equations  of  motion  as 


d0  _ cos(d>) _ 

dh  "  (1  ♦  h)  tan(y)  cos(4>) 


(3.26) 


d$  sin(4>) 
dh  "  (1  ♦  h)  tan(y) 


(3.27) 


dV  _  CD  r,  p,e'n/<  VS  g,  rt 
dh  “  '  2m  sin(y)  "  v(l*hF 


(3.28) 
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2m  sin(y) 


(3.29) 


dy 

dh  ' 


;os(a)  g,  r%  \ 

~  ( 1  +  h)2  V2  tan(y)  +  0  ♦  h)  tan(y) 


d<j)  _  CL  r,  pge'h/<  S  sin(a)  cos(q> )  tan(<)>) 
dh  ’  2m  sin(y)  cos(y)  (1  ♦  h)  tan(y) 

Thus,  the  above  modified  equations  of  motion  are  expressed  in  a 
convenient  integration  parameter,  h.  In  the  next  section,  the  dependent 
variables  are  transformed  into  a  set  of  orbital  parameters  which  are 
convenient  to  use  in  the  analytical  analysis  performed  in  Section  IV. 


Transformation  to  Orbital  Elements 


Ballistic  Coefficient  and  Other  Flight  Parameters.  InEqs(3.26)- 
(3.30),  there  are  several  lengthy  constant  coefficients  which  pertain  to  the 
flight  vehicle  and  the  planetary  atmosphere  it  is  entering.  By  developing 
non-dimensional  shorthand  expressions  for  some  of  these  terms,  the 
equations  become  easier  to  manipulate.  The  first  non-dimensional  term 
defined  is  called  the  ballistic  coefficient,  B,  and  specifies  physical 
characteristics  of  the  flight  vehicle,  as  well  as  the  planetary  atmosphere  it  is 
entering.  B  is  assumed  constant  throughout  the  entry  trajectory,  specified 
by  initial  conditions  and  defined  as  (Vinh  and  others,  1980:256) 


P|SCq  PtrtSCn€ 
Dm  2mp  *  2m 


(3.31) 


To  help  further  facilitate  the  manipulation  of  the  equations  of  motion 
the  terms  specifying  the  bank  angle,  o,  and  the  lift-to-drag  ratio,  CL/CD,  are 
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combined.  6  and  A  are  assumed  constant  throughout  the  entry  trajectory, 
specified  by  initial  conditions  and  defined  as  (Vinh  and  others,  1980:255) 


A  - 


Ci 

^•cos(o) 


(3.32) 


(3.33) 


Substituting  the  above  three  relations  into  Eqs  (3.26M3.30)  gives 

d0  _ cos(4j) _ 

dh  "  ( 1  +  h)  tan(y)  cos($) 


d<|>  sinpjjj 
dh  “  (1  +  h)  tan(y) 


(3.35) 


dV  _  BVe~h/€  g,r« 
dh  '  "  e  sin(y)  '  V(  1  ♦  h)2 


(3.36) 


dy  BA  e~h/€  g«  r,  t  1 

dh  '  e  sin(y)  ’  (l  +  h)2  V*  tan(y) +  ( 1  +  h)  tan(y) 


dd>  B6  e~h/t  cos(4>)  tan(4>) 
dh  ’  c  sin(y)  cos(y)  ’  ( 1  ♦  h)  tan(y) 


Modified  Speed  Ratio.  Recent  analytical  flight  mechanics  studies 
have  determined  that  transforming  the  velocity  terms,  V  and  dV/dh,  into 
the  modified  Chapman  variable,  u,  places  the  equations  of  motion  into  a 
more  practical  form.  By  definition,  u  is  defined  as  (Vinh  and  others, 
1980:229) 
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(3.39) 


V2  cos2(y)  V2(  1  +  h)  cos2(y) 

Solving  for  V  yields 

cos(y)  V  ( 1  ♦  h) 

As  a  convenience  u  is  also  known  as  the  speed  ratio  since  it  relates 
the  vehicle  s  local  horizontal  velocity  to  the  circular  orbital  velocity  at  the 
vehicle's  current  altitude  (Karasopoulos,  1988:45).  To  replace  V  with  u  in 
Eqs  (3.34M3.38),  Eq  (3.39)  is  differentiated  with  respect  to  h,  keeping  g, 
and  rs  constant. 

du  V2  cos2(y)  2V(1  ♦  h)  cos2(y)  dV 
dh  ’  g,r,  *  g,r,  dh 


2V2(1  +  h)  sin(y)  cos(y)  dy 
gsrs  dh 


(3.41) 


Substituting  Eqs  (3.36)  and  (3.37)  into  the  above  equation  transforms 
the  dV/dh  equation  of  motion  into  du/dh  form  and  replaces  Eq  (3.36) 


du  u  2Bu(l+\  tan(y)) 
dh  “  (l*h)  csin(y)  e 


(3.42) 


Modified  Flight  Path  Anale.  In  Eq  (3.37),  sin(y)  terms  appear  in  the 
denominator.  Most  realistic  entry  trajectories  begin  with  -10*  *  y  <  0*  and 
often  result  in  aerodynamic  skipping,  where  y  switches  from  a  negative  to 
positive  quantity.  During  this  transition  from  positive  to  negative  values,  y 
-  0*.  which  could  create  a  singularity  in  Eq  (3.37)  since  l/sin(y),  where  y  - 
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0*.  is  not  defined.  To  eliminate  this  potential  singularity  as  well  as  to 
simplify  the  equations  of  motion,  the  flight  path  angle,  y ,  is  transformed  to 
the  variable  q,  given  as  (Vinh  and  others,  1980:257) 

q  -  cos(y)  (3.43) 

Differentiating  the  above  relation  with  respect  to  h  yields 

dh*'sin(y)<ih  (344) 


Using  the  above  equation  and  the  definition  for  the  speed  ratio 
transforms  Eq  (3.37)  into  a  differential  equation  for  q. 


dj  _g_(  A  Bxe~h/t 
dh  1+hi1  uj  e 


(3.45) 


Thus,  the  equations  of  motion  are  now  given  as 

du  u  2Bu(l+\  tan(y))  h/e 
dh""(l+h)"  esin(y)  e 


(3.46) 


<*..  q  f,  d2)  B*e-h/< 
dh  1  ♦  h  (  u  j  c 


(3.47) 


d0  _ cos(tj)) _ 

dh  ’  (1  ♦  h)  tan(y)  cos($) 


(3.48) 


d£  sin((j>) 
dh  "  ( 1  +  h)  tan(y) 


(3.49) 


dd>  B6  e'h/<  cos(d;)  tan(<t>) 
dh  *  e  sin(y)  cos(y)  ’  ( 1  ♦  h)  tan(y) 


20 


where  q  -  cos(y). 


Classical  Orbital  Elements.  Before  using  the  Method  of  Matched 
Asymptotic  Expansions  to  analytically  solve  the  above  equations,  they  are 
transformed  again  to  express  them  in  a  form  which  will  simplify  the 
analysis  found  in  the  next  section.  Currently,  the  independent  variables  are 
the  speed  ratio  (u),  modified  flight  path  angle  (q),  longitude  (6),  latitude  (4>) 
and  heading  angle  (4)).  During  the  orbital  lifetime  of  a  satellite,  these 
variables  are  constantly  changing.  By  transforming  the  variable  set  (0.  4> 
and  i|> )  to  the  set  of  classical  orbital  elements  (I,  Q  and  a),  the  resulting 
analysis  is  greatly  simplified.  This  simplification  results  from  the  classical 
orbital  elements  being  constant  for  two-body  motion,  where  there  are  no 
perturbing  forces  (Wiesel,  1989:34-35, 58).  As  shown  in  the  next  section, 
the  equations  of  motion  will  be  analyzed  in  two  domains,  exo-atmospheric 
and  atmospheric  flight.  This  transformation  will  greatly  simplify  the  exo- 
atmospheric  (two-body)  analysis,  since  two  of  the  equations  of  motion  are 
constant. 


Inclination  Angle.  1.  To  derive  the  equation  of  motion 
describing  the  inclination  angle,  the  spherical  trigonometric  relationship 
(A.7)  is  used  and  given  below. 

cos(I)  -  cos(<j>)  cos($)  (3.51) 

Differentiating  the  above  equation  with  respect  to  h  and  solving  for 
dl/dh  gives 

dtf  sin(<fr) cos(4>) d&  cos(4>) sin(4;)d<l) 
dh’  sin(I)  dh  *  sin(I)  dh 


21 


Substituting  in  Eqs  (3.49)  and  (3.50)  for  the  differential  relations 
d4>/dh  and  d^/dh  and  using  the  spherical  trigonometric  equations  (A.5)  and 
(A.ll)  gives 


dl_  B  6  cos(a)  _h/€ 
dh  ’  e  sin(y)cos(y) e 


(3.53) 


Longitude  of  the  Ascending  Node.  ft.  To  derive  the  equation  of 
motion  describing  the  longitude  of  the  ascending  node,  the  spherical 
trigonometric  relationship  (A.  15)  is  used  and  given  below. 

sin(<l> )  -  sin(I)  cos(6  -  ft)  (3.54) 


Differentiating  the  above  equation  with  respect  to  h  and  solving  for 
dft/dh  gives 

dft  cos(d))  dd?  cos(6  -  ft)  cos(I)  dl  d£  . 

dh  *  sin(6  -  ft)  sin(I)  dh  "  sin(0  -  ft)  sin(I)  dh  +  dh 


Substituting  in  Eqs  (3.49),  (3.50)  and  (3.53)  for  the  differential 
relations  d<t>/dh,  d<|>/dh  and  dl/dh  and  using  the  spherical  trigonometric 
equations  (A.4),  (A.  10)  and  (A.l  1)  gives 

dft  B  6  sin(q)  _h/<  .. 

dh  "  c  sin(I)  sin(y)  cos(y) e 


Argument  of  Latitude  at  Enoch,  cl  To  derive  the  equation  of 
motion  describing  the  argument  of  latitude  at  epoch,  the  spherical 
trigonometric  relationship  (A.4)  is  used  and  given  below. 


sin,a)  •  Hw 


(3.57) 
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Differentiating  the  above  equation  with  respect  to  h  and  solving  for 


da/dh  gives 

da  cos(<fe)  d4>  sin(a)  cos(I)  dl  ^ 

dh  “  cos(a)  sin(l)  dh  cos(a)  sin(I)  dh 

Substituting  in  Eqs  (3.49)  and  (3.53)  for  the  differential  relations 
d<t>/dh  and  dl/dh  and  using  the  spherical  trigonometric  equation  (A.l  1) 
gives 


da  1  B  6  sin(a)  _h/€ 

dh  "  (1+h)  tan(y) '  e  tan(I)  sin(y)  cos(y) e 


Summary 


The  above  derivations  develop  the  three-dimensional  equations  of 
motion  for  atmospheric  entry.  They  relate  the  classical  orbital  elements  to 
the  non-dimensional  altitude. 


du  _u_  2Bu(l+A  tan(y))  _h/< 
dh  (1+h)  ~  csin(y)  e 


da  _9_fai.  B\e~h/< 

dh  " (1+h) ^  u  j ’  e 

<H  B  6  cos(a)  _h/€ 
dh  "  e  sin(y)cos(y) e 

dfl  B  6  sin(a)  _h 
dh  "  e  sin(I)  sin(y)  cos(y) e 


(3.60) 


(3.61) 


(3.62) 

(3.63) 
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ar|Q 


(3.64) 


1  B  6  sin(ct) 

( 1+h)  tan(y)  e  tan(I)  sin(y)  cos(y) 


,-h/e 
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IV.  Solutions  to  the  Equations  of  Motion  Using  Matched 

Asymptotic  Expansions 


In  Section  HI,  the  five  equations  of  motion  for  three-dimensional 
atmospheric  entry  were  derived  and  are  repeated  below. 


du 
dh  “ 

u  2Bu(l+A  tan(y))  h/- 
"(l+h)~  esin(y)  e 

(4.1) 

$- 

(4.2) 

dh 

(l+h)(u  J  e 

dl 

B  6  cos(a)  _h/. 

(4.3) 

dh  " 

€  sin(y)cos(y) e 

dQ 

B  6  sin(a)  h/. 

(4.4) 

dh  * 

€  sind)  sin(y)  cos(y) e 

da 

1  B  6  sin(a)  h/. 

(4.5) 

dh  * 

(l*h)  tan(y) "  e  tan(I)  sin(y)  cos(y) e 

The  above  five  differential  equations  are  first  order,  non-linear 
ordinary  differential  equations  (ODEs).  Although  there  exist  techniques  to 
solve  any  first  order,  linear  ODE.  the  non-linearity  and  coupling  among  the 
five  equations  prevents  them  from  being  solved  by  traditional  analytical 
means  (Rainville  and  Bedient,  1981:36).  To  preclude  a  numerical  solution 
to  the  above  ODEs,  and  thus  retaining  some  insight  into  the  mechanics  of 
the  problem,  a  higher  order  analytical  solution  is  used  (Vinh  and  others, 
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1980:254).  The  solution  method  implemented  is  called  the  Method  of 
Matched  Asymptotic  Expansions. 

The  Method  of  Matched  Asymptotic  Expansions 

The  problem  of  atmospheric  entry  is  fundamentally  a  situation  where 
an  initially  small  perturbation,  the  atmosphere,  is  introduced  into  the 
equations  of  motion,  but  its  exponentially  increasing  effect  causes  the 
orbital  elements  to  undergo  very  rapid  changes  over  a  very  narrow  region 
of  the  independent  variable,  the  vehicles's  attitude  (Nayfeh,  1981:270-279). 
Thus,  the  atmospheric  entry  equations  are  modeled  as  if  the  atmosphere 
acts  as  a  boundary  layer  adjacent  to  the  orbital  region  solely  specified  by 
Kepterian,  two-body,  motion,  neglecting  other  perturbing  forces  (Vinh  and 
others.  1980:259).  From  the  perspective  of  the  entry  vehicle,  its  orbit  is 
initially  prescribed  by  gravitational  forces,  but  as  its  altitude  decreases,  it 
enters  a  boundary  layer  region  formed  by  the  atmosphere.  The 
aerodynamic  forces  of  lift  and  drag  will  vary  drastically  over  a  small  range 
of  altitude,  as  compared  to  the  mean  orbit  altitude,  and  thus  dramatically 
alter  the  orbital  parameters. 

The  Method  of  Matched  Asymptotic  Expansions  is  chosen  over  other 
analytical  techniques  such  as  the  methods  of  multiple  scales  and  straight¬ 
forward  expansions  since  the  Method  of  Matched  Asymptotic  Expansions  is 
more  adept  at  handling  non-linear  differential  equations  (Nayfeh, 

1981:279).  In  this  method,  the  solution  to  the  problem  is  represented  by 
two  expansions,  each  of  which  is  valid  in  part  of  the  problem  domain,  either 
exo-atmospheric  or  atmospheric  altitude.  Since  the  two  expansions  have 
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some  overlap,  they  are  matched  to  create  one  composite  expansion  valid 
over  the  entire  problem  domain. 

Deriving  the  Outer  and  Inner  Expansions  Using  Mathematica™ 

The  process  of  deriving  the  outer  and  inner  expansions  using 
Matched  Asymptotic  Expansions  is  a  very  laborious  and  tedious  exercise  in 
algebraic  bookkeeping.  To  expedite  these  derivations,  a  computerized, 
symbolic  manipulator,  Mathematica™,  is  used  to  perform  the  outer  and  inner 
expansions.  Using  Mathematical  decreases  the  time  performing  algebraic 
manipulations.  Thus,  more  time  is  spent  in  analysis  of  the  derived  solutions. 
The  computer  code  « J  in  this  study  is  presented  in  Appendix  B,  as  is  the 
methodology  used  «.o  develop  the  code. 

In  the  following  expansions,  two  of  the  ten  expansions  are  performed 
to  illustj  ate  the  techniques  involved.  The  remaining  expansions  are 
presented  without  derivation  and  are  as  given  from  Mathematical  output. 

Outer  Expansions 

The  solutions  developed  for  the  exo-atmospheric  (non-boundary 
layer)  portion  of  the  domain  are  called  the  outer  solutions.  These  solutions 
are  developed  from  asymptotic  expansions  of  the  equations  of  motion  using 
the  small  parameter,  e.  The  outer  solution  variables  are  denoted  by  the 
superscript  M°”  and  are  assumed  as  follows  (Vinh  and  others,  1980:259): 

u°  -  UqOi)  ♦  u,(h)e  ♦  u2(h)e2  ♦  OU3) 
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q°  -  q0(h)  +  qt(h)e  +  q2(h)e2  +  0(e3) 

1°  -  I0(h)  +  Ij(h)e  +  I2(fi)c2  ♦  0(£3) 

Q°  -  Ctyh)  ♦  fi,(h)e  +  02(h)€2  ♦  (Xe3) 
a0  -  c^Ui)  ♦  a,(h)c  +  a2(h)c2  +  0(e3) 

y°  -  y0(h)  +  y,(h)e  +  y2(h)c2  +  0(e3)  (4.6) 

The  approximate  solutions  for  lifting  atmospheric  entry  are  of  order 
€°  (Karasopoulos,  1988:82).  The  solutions  derived  to  c1  and  higher  orders  of 
e  act  as  correcting  factors  to  the  zero  order  solutions  and  account  for 
decreasingly  significant  physical  characteristics  of  the  problem.  This  is 
evident  in  Eq  (4.6),  since  the  solution  order  (u,,  u2,..„  for  example)  is 

multiplied  by  corresponding  powers  of  e.  Thus,  the  zero  order  solution 
gives  the  primary  behavior  of  the  variable  and  higher  order  solutions  add 
corrections  which  progressively  bring  this  approximate  solution  in  line  with 
the  actual  solution. 

du/dh  Outer  Expansion.  This  expansion  is  performed  to  illustrate  the 
techniques  used  in  deriving  an  expansion.  The  original  ODE  is  given  in  Eq 
(4.  Das 

du  _u_  2Bu(l+A  tan(y))  _h/< 
dh*’(l+h)"  esin(y)  e 

Next,  straightforward  expansions  for  the  variables  in  the  above  ODE 
are  defined.  The  expansions  for  u  and  y  are  assumed  in  Eq  (4.6)  while 
templates  for  tan(y)  and  l/sin(y)  are  derived  in  Appendix  C  and  given  as 
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Eqs  (C.l  1 )  and  (C.13).  Substituting  the  expansion  for  y  in  the  two  templates 


gives 


[sin(y0  +  y,€  +  y2e2  +  O(e3))] 

1  yicos(y0)  fyi2( 1  *  C0S^Yq))  y2cos(y0))  2  3 

sin2(y0)  €  2sin3(y0)  sin2(y0)  €  +  0  6 

tan(y0  +  y,€  ♦  y2e2  +  0(e3))  -  tan(y0)  +  — jj— -  e 

cos  (y0) 


cos  (y0)  cos*(y0)j 


Substituting  the  above  expansions  in  the  original  ODE  gives 

du0  du,  du2  ?  u0  ♦  utc  *  u2c2  ♦  0(e3)  2Be‘h/c 
dh  +  dh  €+  dh  €  ”  (Uh)  '  e 


(u0  +  u,€+  u2c2  +  0(c3))[l+K  tan(y0  +  y,e  + y2e2  +  0(e3))l 


x  1+Mtan(y0) 


cos  (y0)  cos  (y0) 


s»n(y0)  sin2(y0) 


x  ^  ♦  0(e3) 


At  this  point,  ail  the  terms  on  the  right  hand  side  are  products  of 

terms  which  are  sums  of  constant  coefficients  and  powers  of  6.  Thus, 

multiplying  out  the  above  equation  and  equating  terms  of  identical  powers 

of  €,  will  give  the  desired  expansions  (ODEs)  of  order  e.  Since  exp(-l/e)  is 

*-h/« 


smaller  than  any  power  of  e  as  e->0, 


0  (Nayfeh,  1981:260).  This 


greatly  simplifies  the  expansion  above  and  is  frequently  used  in  this  study. 
Thus,  the  above  expansion  is  multiplied  out  and  the  order  €  terms  are  as 


follows 


,0 ,  ^ . .  JiL 

c  terms.  dh  (1,h) 


(4.7) 


e  terms: 


dh  (1+h) 


(4.8) 


e  terms: 


dh  "’(Mi) 


(4.9) 


da/dh  Outer  Expansion.  As  derived  by  using  Mathematical  the 
dq/dh  expansions  are  grouped  by  e  order  as  follows: 


e  terms 


.  _  _9b_  fSai  , 

Jk.  "  I  .  L.  ..  ”  » 


dh  1  ♦  h  u, 


lterms.  *Lm_Sa_(*MLm!£*L 1  JUJ^ 
er  ■  dh  i  +  h  u0  \rr  i  +  h  u0  1 


2  2, 


(4.10) 


(4.11) 


terms:  V >,2) 

dn  V  u0  uo  Uft£  w  (u0J  u0‘ 
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dl/dh  Outer  Eioansion.  As  derived  by  using  Mathematics™,  the 
dl/dh  expansions  are  grouped  by  e  order  as  follows: 

o  dIo 

e  terms:  ^  -  0 

l  dIi 

e  terms:  -  0 


2 

€  terms:  ^-0 


(4.13) 

(4.14) 

(4.15) 


dfl/dh  Outer  Expansion.  As  derived  by  using  Mathematics™,  the 
dft/dh  expansions  are  grouped  by  e  order  as  follows: 


o  dfio 

€  terms:  -  0 

i  dQ. 

e  terms:  -  0 


(4.16) 

(4.17) 


2 

€  terms: 


(4.18) 


dq/dh  Outer  Expansion.  As  derived  by  using  Mathematica™,  the 
da/dh  expansions  are  grouped  by  e  order  as  follows: 


e°  terms: 


dap  1 

dh  ’  (1  ♦  h)  tan(y0) 


e1  terms: 


dqi  _  _  y> 
dh  ( 1  ♦  h)  sin2(y0) 


e2  terms: 


dq2  _  yi2cos(y0) _ y2 

dh  ( 1  ♦  h)  sin3(y0)  (1  +  h)  sin2(y0) 


(4.19) 

(4.20) 

(4.21) 


31 


Outer  Eipansion  Solutions 


Terms.  The  complete  set  of  e°  term  outer  eipansion  differential 
equations  are  derived  above  and  are  repeated  below. 


du0  u0 
dh  (1+h) 


(4.7) 


dqp  Qp  [qp2 
dh  "  1  +  h  [  u0 


dl* 

dh 


0 


Bl 

dh 


0 


dap  1 

dh  "  (1  +  h)  tan(y0) 


(4.10) 

(4.13) 

(4.16) 

(4.19) 


Solutions  to  this  set  of  differential  equations  are  derived  in  Appendix 
0  and  are  given  below. 


uo 


dp 


(4.22) 

(4.23) 


Ip  “  C03 


(4.24) 


(4.25) 

(4.26) 
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e1  Terms.  The  complete  set  of  e1  term  outer  expansion  differential 


equations  are  derived  above  and  are  repeated  below. 

du}  ut 
dh  '  “(l*h) 


*Ll  Oo  f2V)i  .  _5i_  (Ssl  . 

dh  "  1  ♦  h  u0  u  2  1  ♦  h  u0 


(4.8) 


(4.11) 


(4.14) 


(4.17) 


dh  (1  ♦  h)  sin2(y0) 


(4.20) 


Solutions  to  this  set  of  differential  equations  are  derived  in  Appendix 
D  and  are  given  below. 


ui'TTh 


(4.27) 


Cni(l*h) 


Cl2Vl  +  h 


or - T 

(2  C0IC02(1  >h))2 


(4.28) 


I|  -  C,3 


(4.29) 


Q,  -  CH 


(4.30) 
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(4.31) 


(c j 2 +  CnC02VCoi  )(l+h)  - 


(C 


ii_ 


00,0,21 


( 1  -  C0,2C02)V-Coi  ♦  2(t*h)  -  C0IC02(l+h)2  +  15 


e2  Terms.  The  complete  set  of  e2  term  outer  expansion  differential 
equations  are  derived  above  and  are  repeated  below. 


du2  u2 
dh  "  ‘(l*h) 


(4.9) 


dq2 

dh  ‘q2 


3q0q|2 

uo 


3q02qiUt 


+ 


0 


-0 


da2_  y,2cos(y0)  y2 

dh  ( 1  +  h)  sin3(y0)  ( 1  +  h)  sin2(y0) 


(4.12) 

(4.15) 

(4.18) 

(4.21) 


Inner  Expansions 

The  solutions  developed  for  the  atmospheric  (boundary  layer)  portion 
of  the  domain  are  called  the  inner  solutions.  These  solutions  are  developed 
from  asymptotic  expansions  of  the  equations  of  motion  using  the  small 
parameter,  e.  The  inner  solution  variables  are  denoted  by  the  superscript. 

"i"  and  are  assumed  as  follows  (Vinh  and  others,  1980:260): 

u1  -  u0(0  ♦  u,(^)e  ♦  u205)c2  ♦  0(e3) 
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q*  -  q0(5)  *  q|(0«  *  q2(?)e2  *  0(e3) 

I‘  - 1„<?)  *  I,(5)e  *  I2(?)«2  *  CKe3) 
d  ■  (ty?)  *  0,(5)e  *  02(5k2  ♦  CKe3) 
a1  -  <^(5)  *  a,(5)e  *  a2(5)e2  ♦  CKe3) 

V  -  y0(5)  *  y,(5)e  *  y2(?)e2  *  CKe3)  (4.32) 

To  derive  the  inner  expansions,  a  new  independent  variable  is 
required  to  force  the  equations  of  motion  to  focus  on  the  boundary  layer 
(Nayfeh,  1981:262).  The  magnified  variable,  becomes  the  new 
independent  variable  for  the  five  coupled  equations  of  motion.  Thus  the 
equations  of  motion  undergo  a  stretching  transformation  which  focuses 
them  on  the  behavior  found  in  the  boundary  layer.  £ .  the  magnified  non- 
dimensional  altitude,  is  defined  as 

S  “  or  h  - c ^ 

d 

Using  the  chain  rule  for  differentiation  gives  an  expression  for  ^ . 

A  dh  d  d  d  X  A  (a 

d^“d?dh“€dh  dh  "  e  d$  (4,33) 


Substituting  the  above  transformations  for  both  h  and  dh  into  the 
original  ODEs  (Eqs  (4.1 )  -  (4.5))  gives  the  ODEs  as  functions  of  the 
stretched  variable,  £ 


du  eu  2Bu(  1  ♦  A  tan(y))  ^ 

dl;  "  "(l+€^)'  sin(y)  e 


(4.34) 
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da 


eg  /q2 
(1  ♦  eJ;)  (u 


-BKe* 


dl  B6cos(g) 

"  sin(y)cos(y) e 

dQ  B6sin(a) 

d^  *  sin(I)  sin(y)  cos(y) e 

da  _ e _  B  6sin(g) 

dl^  ’  (1  +  e^)  tan(y) "  tan(I)  sin(y)  cos(y) e 


(4.35) 

(4.36) 

(4.37) 

(4.38) 


This  expansion  is  performed  to  illustrate  the 


techniques  used  in  deriving  an  expansion.  The  original  ODE  is  given  in  Eq 


(4.34)  as 


du  eu  2Bu(l  ♦  Atan(y))  , 

dJ;  “  "(1+e^)"  sin(y)  e 


Next,  straightforward  expansions  for  the  variables  in  the  above  ODE 
are  defined.  The  expansions  for  u  and  y  are  assumed  in  Eq  (4.32)  while 
templates  for  tan(y),  l/sin(y)  and  1/(1  ♦  e  are  derived  in  Appendix  C  and 
given  as  Eqs  (C.13).  (C.l  1)  and  (C.16),  respectively.  Substituting  the 
expansion  for  y  in  the  two  trigonometric  templates  gives 

[sin(y0  ♦  y,c  ♦  y2e2  ♦  OU3))]-1 

1  y,cos(y0)  fy,2(l  ♦cos2(y0))  y2cos(y0))  2 

■  “ — / — r  -  — — 5 —  e  ♦  - 5 - -  ■  \ —  e  +  0(e  ) 

sin(y0)  sin2(y0)  2sin3(y0)  sin2(y0) 
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tan(  y0  ♦  y^  +  y2e2  ♦  OU3) )  -  tan(y0)  ♦  — ^  “  c 


JyL«yfi)+  *_U 

cos3(y0)  cos2(y0)J 


Substituting  in  the  expression  (1  ♦  eO  into  the  algebraic  template  gives 
j7^-I-5e*?2c2.0(e3) 


Substituting  the  above  expansions  in  the  original  ODE  gives 

du,  du,  du*  ,  u0*u,e>u2e2.Q(e3) 

d^  d^  6  dJ;  €  e  1  "^e 


(u0*  u,e+  u2c2  +  0(c3))[|*K  tan(y0*  y,e*  y2€2  +  0(e3) ) 
sin(y0  ♦  y,e+  y^*  CKe3)) 


+  0(e3) 


-  -  €  (u0  ♦  u,c  ♦  u2c2  +  0(e3)  )(l-^e  +  ^2e2  +  0(c3))  -  2B  e'^ 
x  (u0  +  u,e+  u2c2  ♦  0(e3)) 


,  y i  fyizsin(y0)  y2  2 

I  "  Yo  +  cos2(y0)  €+  cos  (y0)  cos^j) 6 


I  yicos(y0)  Y|2(lW(y0))  y2cos(v0)  ,  3 

Lsin(yo)  sin2(y0)  2sin3(y0)  sin^(y0)  J 


At  this  point,  all  the  terms  on  the  right  hand  side  are  products  of 
terms  which  are  sums  of  constant  coefficients  and  powers  of  e.  Thus, 
multiplying  out  the  above  equation  and  equating  terms  of  identical  powers 
of  e,  will  give  the  desired  expansions  (ODEs)  of  order  €. 


€°  terms: 


du0  2Bu0e~*(l  ♦  \tan(y0)) 
d^  " '  sin(y0) 


(4.39) 


,  du,  fui(l  +  Atan(y0)) 

e  terms:  -  -u0  -  2Be^  - - 


u0ylcos(y0)(l  >Atan(y0))  A  u0  y, 

sin2(y0)  cos2(y0)  sin(y0) 


(4.40) 


e2  terms:  ^T?*-u,*u0?- 2Be-<  [-£- * 

d?  0  [sin(y0)  ^cosz(y0)  cosJ(y0) 


y,  ^ 1  ♦  cos  l 
u0(.  ^wn(y0)|— — 3-- 


2(l  ♦  cos2(y0)) 


y2cos(y0) 

sin2(y0) 


Auiy,  ulyicos(y0)(l  +  Atan(y0)) 
cos2(y0)  sin(y0)  sin2(y0) 


^UpYi 

sin2(y0)  cos(y0) 


u2(l  ♦  Atan(y0)) 
sin(y0) 


(4.41) 


do/dE  Inner  Expansion.  As  derived  by  using  Mathematical  the 
dq/dh  expansions  are  grouped  by  c  order  as  follows: 
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e°  terms:  -  -  \  Be"^ 


(4.42) 


e1  terms:  -5^  -  q0 


-1  ♦ 


V 


a? 

Ua 


(4.43) 


dq2  q03  /  \ 

terms:  ~rf  •  (l-u,^u0J  ♦  q0(^-l) ♦  q, 

US  Un 


2  \ 


-  1 


Hfll 


luo  ) 


2qpQtQ2 


u. 


(4.44) 


dl/d l  Inner  Expansion.  As  derived  by  using  Mathematical  the  dl/dh 
expansions  are  grouped  by  e  order  as  follows: 


ft  dl0  B6e~^  cos(oa) 

e  terms:  df-sin(y0)cos(^) 


(4.45) 


dl 


e1  terms:  -  B6e 


-f. 


yjCOsCc^)  (• 


1 


1 


cos2(y0)  sin2(y0)j 


aiSinCop) 
sin(y0)  cos(y( 


ol 


e2  terms:  ^-B6e  *  (^costc^) 


(4.46) 


yi2cos(o<o)  /  1  1 

sin(y0)  cos(y0)  (cos2(y0)  sin2(y0) 

(~2“  ♦  y,2]cos(a0)  ♦  otj  sinCap) . 
sin(y0)  cos(y0) 


(4.47) 
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dfl/dg  Inner  Expansion.  As  derived  by  using  Mathematical  the 
dQ/dh  expansions  are  grouped  by  e  order  as  follows: 


e  terms: 


dflp  B6e~*  sinfcip) 


d£  sin(I0)  sin(y0)  cos(y0) 


(4.48) 


i  _ ryisin(ao)  /  i 


y0)  sin2(y0) 


If  sinCap)  cos(IQ)  oncost  otp) 

sin2(I0)  sin(y0)  cos(y0)  sin(I0)  sin(y0)  cos(y0) 


(4.49) 


e2  terms:  ^  -  B6e‘*  [ ( — t. — 7  -  — r: — I" 


cos2(y0)  sin2(y0)JL  sin(I0) 


y.I.sintoo)  cos(I0) 


sin(y0)  cos(y0)  sin(ao)  sin(I0) +  sinCoio)  sin(I0)  ‘  sin2(I0) 


1 

fi,li 

(l  +cos2(I0)) 

sin(y0)  cos(y0)  sinfctp)  sin(I0) 

2sin2(I0) 

cos(IQ) 

sin(I0) 


(l2  ♦  a^sinfap)  cos(ao))J| 


(4.50) 


da/dg  Inner  Eioansion.  As  derived  by  using  Mathematical  the 
da/dh  expansions  are  grouped  by  e  order  as  follows: 
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0  do,  B6e~*  sinCop) 
e  terms.  ^  sjn(yo)  Cos(y0) 


(4.51) 


e1  terms: 


det|  1  nr 

■y,sin(op) 

(  1  1  \ 

d(|  tan(y0)  D  * 

tan(I0) 

cos2(y0)  sin2(y0); 

I,  sin(op) 

a,cos(op) 

(4.52) 


e  terms: 


da 


Yi 


L 


^  sin2(y0)  tan(y0) 


-  B6e 


(cosz(y0)  sinz(y0)JL  tan(I0) 


y1a1cos(a0) 


_ ll _ + _ V2 

sin(y0)  cos(y0)  sintap)  tan(I0)  sin(op)  tan(I0) 


y^sintop) 

sinZ(I0) 


+ 


1 

I,2  cos(Ip) 

sin(y0)  cos(y0)  sin(ap)  tan(I0) 

.  2sinz(I0) 

+  sinz(a)) 


cos(I0)cos(q0)(a2cos(a0)-q12sin(a0))  -  ^  qlIls^VSPs(qQ) 


(4.53) 


Inner  Expansion  Solutions 

e°  Terms.  The  complete  set  of  e°  term  inner  expansion  differential 
equations  are  derived  above  and  are  repeated  below. 

du0  2Bu0e~*(l  ♦  htan(y0)) 
d^  “ '  sin(y0) 
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(4.39) 


dq0 


-\Be~* 


dl0  B6e~*  005(0(0) 

"  sin(y0)  cos(y0) 

dflp  B6e~*  sin(qp) 
dJ;  “  sin(I0)  sin(y0)  cos(y0) 

dotp  B6e~*  sin(qp) 

d^  “  ’  tan(I0)  sin(y0)  cos(y0) 


(4.42) 

(4.45) 

(4.48) 

(4.51) 


Solutions  to  this  set  of  differential  equations  are  derived  in  Appendix 
E  and  are  given  below. 

uo  “  Koi(*Be'*  ♦  K02)2  expj^Cos'1  (\Be“*  ♦  K02)j  (4.54) 


q0  -  ABe'?  ♦  K02 
I0  -  Cos*1  cos(K03)cos[^  ln|tan  j  ♦ 


'ff  Cos',(\Be”^  +  K02)l 


+  K, 


05 


(4.55) 

(4.56) 

(4.57) 


-  Cos 


-1 


H 


sin2(K03) 


1T\  %\ 


Ll-cos2(IC03)cos 


n  Cos"1  (\Be~?+K 


02, 


Me. 


05 


V 


cos(K03) 


} 


dp  -  Sin"1 


/sindC03)\ 

( sin(I0)  J 


(4.58) 
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-  Sin' 


sinOCflJ 


i  ?ri'  \  2  6  i  [,  f77  C°S  ^Be  5+IC°2^  ))  1]/2 

1-cos4(K03)  cos|^  in^tan^+ - 5 - J]+K05 


ej_  Terms.  The  complete  set  of  e1  term  inner  expansion  differential 
equations  are  derived  above  and  are  repeated  below. 

du,  Fu,(l  +  Man(y0)) 

"  U°  "  2Be  sin(y0) 


u0yic°s(y0)( 1  +  ^fanfyo))  \  u0  y. 

- 5 - +  — 5 — ^ - 

sin  (y0)  cos  (y0)  sin(y0). 


(4.40) 


*!l  f,  3 l 

-‘‘t- 


dl, 

df'B5e 


5  '  sin2(y0)) '  sin(y0)  cXoJ 


(4.43) 


(4.46) 


§l.B6e-JVl^(  ■  .  1_ 

d^  [  sm(I0)  (cos2^)  sin2(y0) 


I.  sin(oft)  cos(Ift) 


a,cos(cu) 


sin2^)  sin(y0)  cos(y0)  sin(I0)  sin(y0)  cos(y0) 


(4.49) 
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d£l.  _1 _  -)[V'ii"Wf  1  1  \ 

d?  tan(y0)  land,))  |cos2(y0)  sin2(y0)j 


It  sinCop)  o^costap) 

sin2(I0)  sin(y0)  cos(y0) +  kn(I0)  sin^o)  C0S<V 


(4.52) 


Solutions  to  this  set  of  differential  equations  are  derived  in  Appendix 
E  and  are  given  below. 

u,  -  (ABe‘?*  K02)2eipf-^Cos'‘(\Be"?*  K02)| 

0  .  2KI2I  K0I  Sin-'fABe-^K  ) 
l  "  AB<f<  *  K02  A  Sl"  'ABe  Ko2> 

h  1*01^12 ~  *01  *ml(*«25-ABe  Q  2K0;  K„|  K|2i 
\jl .  (^Be^-Koj)2  W>  -Ka22 

j  2^(1-K022)  [  l-(ABe~;»K02)2 ]  -  2K02(ABe~^  K02)  ♦  2  . 

L  ABe'2  JJ 


(4.59) 


i.-k,2,  (k02<;-ab«-')*k12 


(4.60) 


44 


45 


2  7 

f_-Terms.  The  complete  set  of  c  term  inner  expansion  differential 
equations  are  derived  above  and  are  repeated  below. 


^  s,n^o)  (cos2(y0)  cos3(y0) 


u0(l*  ^n(Y0))[Y'^1*3Os2(Y0^  - 

01  0  '[  2sin3(y0)  sin2(y0) 


KUiVi  ulylcos(y0)(l  ♦  Ntan(y0)) 


♦  5  " 

cos  (y0)  sin(y0) 


sin2(y0) 


^pyi2  u2(l  +  htan(y0)) 
sin2(y0)cos(y0)+  sin(y0) 


(4.41) 


(4.44) 


B6e  5  (y2cos(oo)  -  yja^infotp)) 


cos2(y0)  sin2(y0) 


y,  cos(ao)  /  i  ] 

sin(y0)cos(y0)  cos2(y0)  sin2(y0)  1 


a.  2 

2  +  Yj  Jcosfop)  ♦  ajsi^ap) 
sin(y0)cos(y0) 


(4.47) 


46 


-  B6e"? 


‘y^cosCap) 

sin(I0) 


_ Yl _  y2  yiIlsin(a0)  cos(I0) 

+  sin(y0)  cos(y0)  510(00)  sin(I0) +  sin(oo)  sin(I0) "  sin2(I0) 


l 

I,2(l  *cos2(I0)) 

sin(y0)  cos(y0)  sinCo^)  sin(I0) 

2sin2(I0) 

a2sin(ao)  cos(oo) 


*  otjl.sinfoto)  costal)  J 


(4.50) 


da 


_lL 


JL 


^  ’  "  sin2(y0)  ’  tan(y0) 


-  B6e 


-t,  If  1  l  xTyia^ostoto) 

,cos2(y0)  sin2(y0)JL  t«n(I0) 


+ _ Yl _  y2  y^sinCotp) 

+  sin(y0)  cos(y0)  sin(ao)  tan(I0) +  sin(oo)  tan(I0) "  Sin2(I0) 


1 

l,2cos(Ip) 

sin(y0)  cos(y0)  sin(oo)  tan(I0) 

.  2sinz(  I0) 

+  cos(l0)cos(a0)(a2cos(a0)-a,2sin(a0))  -  12-qlIl^^!C0SfqQ^  (4.53) 


Matching  Asymptotic  Eioansions 
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Once  the  asymptotic  expansions  for  the  equations  of  motion  are 
derived,  the  outer  and  inner  solutions  are  matched*  to  reduce  the  number  of 
unknowns,  the  constants  of  integration,  resulting  in  an  equal  number  of 
initial  conditions  and  constants  of  integration.  Thus,  the  final  solutions  to 
the  equations  of  motion,  known  as  the  composite  solutions,  are  used  to  solve 
entry  trajectories. 

The  underlying  idea  of  using  outer  and  inner  solutions  is  to  derive 
two  solutions  which  model  the  two  dominant  forces  acting  on  the  entry 
vehicle,  gravity  and  air  drag.  To  model  each  of  these  forces,  two  scales,  the 
independent  variables,  are  used  to  derive  two  separate  expansions;  each  are 
valid  in  part  of  the  domain  of  the  entry  problem,  but  neither  covers  the 
entire  domain  of  interest.  Although  each  of  the  scales  do  not  cover  the 
entire  domain,  they  do  overlap,  or  have  regions  where  both  expansions  are 
valid.  Since  the  expansions  having  neighboring  regions  of  validity,  they  are 
blended  or  matched,  resulting  in  one  composite  solution  which  connects  the 
two  previously  separate  solutions. 

In  this  study.  Van  Dyke's  matching  principle  (Nayfeh,  1981:282-283) 
is  used  to  equate  the  outer  and  inner  constants  of  integration.  Van  Dyke  s 
principle  states  that  the  m^  inner  expansion  of  the  0th  outer  expansion 
equals  the  n^  outer  expansion  of  the  mth  inner  expansion,  where  m,n  are 
integer  values.  Simply  stated,  where  y  is  the  derived  expansion  and  i  and  o 
are  the  outer  and  inner  expansions 

[(y)°]'-[(y)T  «M) 


48 


To  obtain  the  inner  expansion  of  the  outer  expansion,  the  outer 
expansion  is  rewritten  in  terms  of  the  inner  variable,  J; ,  and  is  expanded  for 
small  e ,  keeping  ^  constant.  To  obtain  the  outer  expansion  of  the  inner 
expansion,  the  inner  expansion  is  rewritten  in  terms  of  the  outer  variable,  h. 
and  is  expanded  for  small  e ,  keeping  h  constant.  The  e°  terms  of  the  two 
resulting  expressions  are  equated,  thus  expressing  the  inner  constants  of 
integration,  Kjj ,  in  terms  of  the  outer  constants  of  integration,  Ciy  correct  to 

CKc). 


Zero  Order  c  Matching 

In  the  following  sections,  Van  Dyke's  matching  principle  is  applied  to 
the  c°  inner  and  outer  expansion  solutions.  The  outer  expansions  will  be 
rewritten  in  terms  of  the  inner  variable  first,  and  then  the  inner  expansions 
will  be  rewritten  in  terms  of  the  outer  variable.  Both  composite  expansions 
will  be  expanded  for  e  and  their  resulting  e°  terms  will  be  equated,  resulting 
in  the  inner  expansion  constants  of  integration  being  expressed  in  terms  of 
the  outer  expansion  constants  of  integration. 

Inner  Expansion  of  the  Outer  Expansion  Solutions.  In  this  section,  the 
five  outer  expansion  solutions  are  rewritten  in  terms  of  the  inner  variable 
and  expanded  to  obtain  c°  terms. 

Matching  Expansion  for  u0.  From  Eq  (4.22),  the  c°  outer 
expansion  solution  for  u  is 


Rewriting  the  outer  variable,  h,  in  terms  of  the  inner  variable.  I ,  uses 
the  definition  h  -  e  ^ .  Substituting  this  in  the  above  equation  gives 

/u  \omS&L- 

W  l*e$ 

Using  the  binomial  expansion  derived  as  Eq  (C.16)  gives 
J7J7-  1  -Se.^.OU3)-  1  *0(«) 

Thus,  the  inner  expansion  of  the  e°  outer  expansion  solution  for  u  is 

[(“o)0!* -coi  (4.65) 


From  Eq  (4.23),  the  eu  outer 


expansion  solution  for  q  is 


(n  \n  C0lV^02 

0  Vi  -|colCo2(i  *h)-i|: 


2d  ♦  h)  -  Col  C02d  .  h)2 


Rewriting  the  outer  variable,  h,  in  terms  of  the  inner  variable,  i; ,  uses 
the  definition  h  -  e  ^ .  Substituting  this  in  the  above  equation  gives 


VCqjl 


(q  \0- _ _ 


Using  the  binomial  expansion  derived  as  Eq  (C.17)  gives 


(-C0iC02£252 .2d  -  c0lc02)«5 .  (2  -  Co.Coj))*4  "  (2  -  C0IC02)' 
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Thus,  the  inner  eipansion  of  the  e°  outer  expansion  solution  for  q  is 


[Mi- 


'M. 


2  ■  Cn  i  C, 


oi  M)2 


(4.66) 


Matching  Expansion  for  I0.  From  Eq  (4.24),  the  e°  outer 
expansion  solution  for  I  is 

0o)°  "  ^03 

By  inspection,  since  there  is  no  dependence  on  h,  the  inner  expansion 
of  the  e°  outer  expansion  solution  for  I  is 


[Mi- 


'03 


(4.67) 


Matching  Expansion  for  ftp.  From  Eq  (4.25),  the  e°  outer 
expansion  solution  for  Q  is 

(flof-C* 

By  inspection,  since  there  is  no  dependence  on  h,  the  inner  expansion 
of  the  e°  outer  expansion  solution  for  ft  is 


[W]  *C<H 


(4.68) 


Matching  Expansion  for  oq..  From  Eq  (4.26),  the  e°  outer 
expansion  solution  for  a  is 

C. 


(a,,)0  -  Sin*1 


1  - 


1  *h 


IT- 


C0|2C02J 


♦C 


05 
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Rewriting  the  outer  variable,  h,  in  terms  or  the  inner  variable,  1; ,  uses 
the  definition  h  -  e  ^ .  Substituting  this  in  the  above  equation  gives 


(ao)°'Sin  '  ylic*c  tCoi 

L  V  1  ■  '"'01  LD2J 


Using  the  binomial  expansion  derived  as  Eq  (C.16)  gives 


1  +  e  ^ 


-  1  ♦  0(e) 


Thus,  the  inner  expansion  of  the  eu  outer  expansion  solution  for  a  is 


[  (ao)0]' '  Sin' 


‘-Cqi 

'l  -  cJc 


(4.69) 


01  ^02-J 


In  this  section,  the 

five  inner  exoansion  solutions  are  rewritten  in  terms  of  the  outer  variable 
and  expanded  to  obtain  e°  terms. 


j.  From  Eq  (4.54),  the  €°  inner 


expansion  solution  for  u  is 


K)1  -  K01(ABe-« .  K„,)2  eip[-f  Cos-'()iBe-{  *  K0J)j 


Rewriting  the  inner  variable,  ^ ,  in  terms  of  the  outer  variable,  h,  uses 
the  definition  l  -  h/e .  Substituting  this  in  the  above  equation  gives 

(u0)‘  -  K0l(\Be-h/‘  *  K02)2eip[-£cos-|(*Be-h/< .  K^)] 
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Since  exp(-l/e)  decreases  to  zero  rapidly  as  e-»0,  e‘h/<  -  0  (Nayfeh, 
1981:260).  Therefore 

\  Be~h/e  ♦  (C02  ■  ^02 

Thus,  the  outer  expansion  of  the  e°  inner  expansion  solution  for  u  is 

[(u0)‘]°  -  K01K022  expj^Cos‘‘(K02)]  (4.70) 

Matching  Expansion  for  o0.  From  Eq  (4.55),  the  e°  inner 
expansion  solution  for  q  is 

(doV-KB^Koj 

Rewriting  the  inner  variable,  ^ .  in  terms  of  the  outer  variable,  h,  uses 
the  definition  ^  -  h/e .  Substituting  this  in  the  above  equation  gives 

(q0)i-*Be-h/,.K02 

Since  exp(-l/e)  decreases  to  zero  rapidly  as  e->0,  e‘h/<  -  0  (Nayfeh, 
1981:260).  Therefore 

^e'h/<*K02-K02 

Thus,  the  outer  expansion  of  the  e°  inner  expansion  solution  for  q  is 

luvr-Koi  <47» 

Matching  Expansion  for  I0.  From  Eq  (4.56),  the  c°  inner 
expansion  solution  for  I  is 
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Rewriting  the  inner  variable,  ^ ,  in  terms  of  the  outer  variable,  h,  uses 
the  definition  ^  -  h/e .  Substituting  this  in  the  above  equation  gives 
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(Oo)*  •  Km 


Since  exp(-l/e)  decreases  to  zero  rapidly  as  e-*0,  e*h/<  -  0  .  Therefore 
\Be“h/e*K02-K02 

Thus,  the  outer  expansion  of  the  e°  inner  expansion  solution  for  0  is 


now  used  to  correlate  the  inner/outer  expansions  and  outer/inner 
expansions.  From  Eq  (4.64)  the  matching  principle  states,  where  Y  is  an 
arbitrary  expansion 

[wf-[<.>T 

Thus,  the  inner/outer  expansions  and  outer/inner  expansions  derived 
above  are  equated,  resulting  the  the  original  inner  expansion  constants  of 
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integration,  K0i,  being  expressed  in  terms  of  the  original  outer  expansion 
constants  of  integration,  C0i.  These  outer  expansion  constants  of  integration 
are  derived  from  the  initial  conditions  of  the  planetary  entry  problem  being 
examined. 


Blending  a0  Matching  Expansions.  Applying  the  matching 

principle,  Eq  {4.64),  to  the  inner/outer  and  outer/inner  expansions,  Eqs 
(4.66)  and  (4.71),  for  q0  gives 


(4.75) 


Blending  u0  Matching  Expansions.  Applying  the  matching 

principle  to  the  inner/outer  and  outer/inner  expansions.  Eqs  (4.65)  and 
(4.70),  for  u0  gives 

Coi-Ko.K^eipf-fcos-'OCoj)] 


Solving  for  K0,  gives 


Koi  -  <2  -  Q 


0,Co2)  ex 


2 

p  ^Cos_,(K02) 


(4.76) 


Blending  Oq  Matching  Expansions.  Applying  the  matching 

principle  to  the  inner/outer  and  outer/inner  expansions,  Eqs  (4.69)  and 
(4.73),  for  Oq  gives 


Sin 


-1 


1  ■  Cf 


sin(Koa) 


;  I  ffl  —  ♦Cab  -Sin'1  •  (f/r  \i  1°1 

bF^Xj Cos  Hiw  1  }J 
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Using  the  matching  principle  and  Eq  (4.67),  C03  -  [  (Iq)0]*"  [Oo)']° 


Substituting  this  in  the  above  relation  and  solving  for  K03  gives 


1C03  -  Sin'1  sin(C03)  sin  Sin 


.  -i  (  1  ~  Qi  ~ 

7 - T=  +  c05 

ky  1  -  c0,  c02> 


(4.77) 


Applying  the  matching 


principle  to  the  inner/outer  and  outer/inner  expansions,  Eqs  (4.67)  and 
(4.72),  for  I0  gives 


C03  -  Cos"1  cos(K03)cos[|  Indian  ^  + 


Cos_,(KnJ 


Solving  for  K05  gives 


ICak  ■  Cos 


l|~COS(CQ3)]  SJu>tC0S'!(ICQ2) 


*5  ^  cosdCn,)  Kmr,,UT  2 


(4.78) 


Blending  Matching  Expansions.  Applying  the  matching 

principle  to  the  inner/outer  and  outer/inner  expansions,  Eqs  (4.68)  and 
(4.74),  for  gives 


C<)4  "  ^04  ‘  Cos 


cos(K03)  ) 


Using  the  matching  principle  and  Eq  (4.69) 


[(ao)°]'-[(ao),]°-sin'1 


■-Co.  1  c 

I  2  L05 

1  ~  C01  C02. 


Substituting  this  in  the  above  relation  and  solving  for  K0<4  gives 
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♦  Cos"1 


Sin'1 

‘-Co.  1 

r  ft 
+  C05 

✓ 

cosJ 

1 - 

1 

o 

o 

O 

o 

KJ 

1 _ 

cos(K03) 

/ 

(4.79) 


Zero  Order  e  Solutions  to  the  Equations  of  Motion 


The  above  derivations  provide  two  separate,  outer  and  inner, 
expansions;  each  is  valid  in  part  of  the  altitude  domain,  but  neither  is  valid 
over  the  entire  domain.  Additionally,  since  the  altitude  where  the  inner 
expansion  is  more  accurate  than  the  outer  expansion  is  not  precisely  known, 
where  to  switch  from  the  outer,  gravity  dominated,  expansion  to  the  inner, 
air  drag  dominated,  expansion  is  a  not  known.  To  generate  an  expansion  or 
solution  valid  over  the  entire  domain,  the  inner  and  outer  expansions  are 
blended  to  form  a  composite  expansion.  Using  this  composite  expansion 
negates  the  requirement  to  switch  from  the  outer  to  inner  expansions,  at  an 
ambiguous  altitude,  to  obtain  the  solution  to  the  planetary  entry  problem. 
The  composite  expansion  for  e°  solutions  is  defined  as  (Nayfeh,  1981:277) 


y\y' 
Yo  +  ‘o 


(4.80) 


The  outer/inner  (or  inner/outer)  expansion  above  accounts  for  the 
components  common  between  the  inner  and  outer  expansions.  Thus,  in  the 
altitude  domain  controlled  by  the  outer  expansion,  the  inner  expansion  is 
negated  by  the  outer/inner  expansion  and  vice  versa.  The  two  possibilities 
above  are  equivalent  as  defined  by  the  matching  principle,  Eq  (4.64). 
Throughout  the  composite  expansions  below,  [yJ]  is  used,  unless  otherwise 
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noted,  since  it  typically  is  a  simpler  and  more  compact  expression  than 

KT- 


u0  Composite  Expansion.  Applying  the  composite  expansion 

definition,  Eq  (4.80),  to  the  outer,  inner  and  inner/outer  expansions,  Eqs 
(4.22),  (4.54)  and  (4.65)  respectively,  for  u0  gives 


uo  *  'coi]Th  +  Koi  ( *Be~h/e+K02  )2  exp^-j"  Cos'1  ( M5e-h/€+K02 ) 


(4.81) 


where  K0I  and  K02  are  given  by  Eqs  (4.76)  and  (4.75). 


q0  Composite  Expansion.  Applying  the  composite  expansion  definition 
to  the  outer,  inner  and  outer/inner  (here  this  expansion  is  simpler  than  the 
inner/outer  expansion)  expansions,  Eqs  (4.23),  (4.55)  and  (4.71) 
respectively,  for  q0  gives 


VHc 


C0l'^02 


oiC02(l  ♦ 


h)- l|: 


♦  A  Be_h/< 


(4.82) 


10  Composite  Expansion.  Applying  the  composite  expansion  definition, 

to  the  outer,  inner  and  inner/outer  expansions,  Eqs  (4.24),  (4.56)  and  (4.67) 
respectively,  for  I0  gives 


c  I 
I0  -  Cos  1 


cosOC03)cos  r  In  tan 


n  Cos’ 1  (  \Be-h/<  +  K02  y 


♦K, 


05 


(4.83) 


where  K02,  K03  and  K05  are  given  by  Eqs  (4.75),  (4.77)  and  (4.78), 
respectively. 
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Applying  the  composite  expansion  definition 


to  the  outer,  inner  and  inner/outer  expansions,  Eqs  (4.26),  (4.58)  and  (4.69) 


respectively,  for  gives 


where  K02,  K03  and  K05  are  given  by  Eqs  (4.75),  (4.77)  and  (4.78), 
respectively. 


Applying  the  composite  expansion 


definition  to  the  outer,  inner  and  inner/outer  expansions,  Eqs  (4.25),  (4.57) 
and  (4.68)  respectively,  for  Qq  gives 


(4.85) 
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where  K02,  K03,  K04  and  K05  are  given  by  Eqs  (4.75),  (4.77),  (4.79) 
and  (4.78),  respectively. 

First  Order  e  Matching 

In  the  following  sections,  Van  Dyke's  matching  principle  is  applied  to 
the  e1  inner  and  outer  expansion  solutions.  The  outer  expansions  will  be 
rewritten  in  terms  of  the  inner  variable  first,  and  then  the  inner  expansions 
will  be  rewritten  in  terms  of  the  outer  variable.  Both  composite  expansions 
will  be  expanded  for  e  and  their  resulting  e°  terms  will  be  equated,  resulting 
in  the  inner  expansion  constants  of  integration  being  expressed  in  terms  of 
the  outer  expansion  constants  of  integration. 

Inner  Expansion  of  the  Outer  Expansion  Solutions.  In  this  section,  the 
five  outer  expansion  solutions  are  rewritten  in  terms  of  the  inner  variable 
and  expanded  to  obtain  e°  terms. 

Matching  Elpansiwifor  Up  From  Eq  (4.27),  the  e1  outer 
expansion  solution  for  u  is 

Rewriting  the  outer  variable,  h,  in  terms  of  the  inner  variable,  Z, ,  uses 
the  definition  h  -  e  Z, .  Substituting  this  in  the  above  equation  gives 


Using  the  binomial  expansion  derived  as  Eq  (C.16)  gives 
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(4.86) 


Y7T7-  CKe3) -  l*0<e) 

Thus,  the  inner  expansion  of  the  e1  outer  expansion  solution  for  u  is 

l(«i)°]i-c„ 


From  Eq  (4.28),  the  e1  outer 


expansion  solution  for  q  is 


(2 -C01C02(1  -  h))2 


Rewriting  the  outer  variable,  h,  in  terms  of  the  inner  variable,  ^ ,  uses 
the  definition  h  -  e  ^ .  Substituting  this  in  the  above  equation  gives 


(2-C01C02(l*e5))2 


Using  the  binomial  expansion  derived  as  Eqs  (C.9)  and  (C.17)  gives 
/  1  -  1  ♦  0(e)  and  \j\  ♦  el;  -  1  ♦  0(e) 

Vl  ♦ 

3  3 

(2  -  C„,C02(  1  *  e  5))'2  -  ((2  -  C0IC02)  -  C0IC02  e  ?)‘2 

00 

■  Z  ('3^)(2-C0|C02)'^r‘(-C0|C02«^ 

k*0 
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00 


[-if  P/A  (2-C0iC02)  3'2'k 


(-CoiCoae^ 


■  (2  -  c01c02)  0(e) 


Thus,  the  inner  expansion  of  the  e1  outer  expansion  solution  for  q  is 


(2  -  C0IC02)2 


(4.87) 


expansion  solution  for  I  is 


(I|)°-C|, 


L,.  From  Eq  (4.29),  the  e1  outer 


By  inspection,  since  there  is  no  dependence  on  h,  the  inner  expansion 


of  the  e1  outer  expansion  solution  for  I  is 


(o.rj'-c,, 


(4.88) 


lv  From  Eq  (4.30),  the  e1  outer 


expansion  solution  for  ft  is 


By  inspection,  since  there  is  no  dependence  on  h,  the  inner  expansion 


of  the  e1  outer  expansion  solution  for  ft  is 


[(ft,)0]  -CM 


(4.89) 
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Rewriting  the  outer  variable,  h,  in  terms  of  the  inner  variable,  Z, .  uses 
the  definition  h  -  e  ^ .  Substituting  this  in  the  above  equation  gives 
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outer  hx pan sion  of  the  inner  Expansion  Solutions.  In  this  section,  the 
five  inner  expansion  solutions  are  rewritten  in  terms  of  the  outer  variable 
and  expanded  to  obtain  e°  terms. 


L,.  From  Eq  (4.59),  the  e1  inner 


expansion  solution  for  u  is 


(u,)1  -  (\Bc  *  +  K02)2  exp|- ^Cos'^Be-^  +  K02)j 
x  (Ktl  >  ^  k02  k12i)  (<  ^  j2  2KiZl  ICgi  ./  g  \ 


K  [K0IKI2  •  <01  Ki2il(K02^ '  2Kn,  Kft,  K, 


l-(\Be-^K02)2 


XVI  -K, 


i 


2M(l-K022)[HNBe‘^ICo2)2l  -  2K02(\Be~*  -  K02)  ♦  2 


XBe‘< 


Rewriting  the  inner  variable,  ^ ,  in  terms  of  the  outer  variable,  h,  uses 
the  definition  ^  -  h/c .  Substituting  this  in  the  above  equation  gives 

(u| )‘  ■  ( \Be-h/<  ♦  IC02)2  exp|-  ^ Cos' 1  ( \Be‘h/t  ♦  K02)j 

X  [k,,  ^  2tC0l(K|2  ^  K02  K12l)  ^  *  l2  ./  h/e  \ 

XBe"h/<  +  Kf  '  X  Sm  *** 


\  [<0I<I2  "  <01  <121  |  (<02  $  '  ^Be  h/<) 
"\j 1  ‘  (^>Be"h/' »  K02)2 


2Kq2  KgjJCj 
M  1  *  Krti 
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2-\j 

f<I-Ko22> 

l-(\Be~h/MC02)2 

■  2K02< 

(\Be-h/,.K0j) 

1  *  2 

1 

KB 

J 

Since  exp(-l/c)  decreases  to  zero  rapidly  as  c-*0,  e*h/e  -  0  (Nayfeh, 
1981:260).  Therefore 

^Bc'h/<+  K02  -  K02 


Thus,  the  outer  eipansion  of  the  e1  inner  expansion  solution  for  u  is 


[(ul)’J°  -  K022  «P 


1  r  Cos‘i(K02)] 


K*l  +  2K«»fl^+IC>2')<5*  l> 

Vn02  / 


*  "  02  ^hTa, 


^J2  S  *  ^121 


M(l-K022)^ 


KB 


(4.91) 


Matching  Expansion  for  ot.  From  Eq  (4.60),  the  e1  inner 
expansion  solution  for  q  is 


Rewriting  the  inner  variable,  ^ ,  in  terms  of  the  outer  variable,  h,  uses 
the  definition  1;  -  h/e .  Substituting  this  in  the  above  equation  gives 

(lt)i-K|2i(K,2?-''Be-|,/')*Ku 

Since  exp(-l/c)  decreases  to  zero  rapidly  as  e->0,  e*h/<  -  0  (Nayfeh, 
1981:260).  Therefore 

KBe"h/<-0 

Thus,  the  outer  expansion  of  the  c1  inner  expansion  solution  for  q  is 
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|(<l|)i]0-K,21Ko2^K12 


(4.92) 


Matching  Eioansion  for  From  Eq  (4.61),  the  e1  inner 
expansion  solution  for  I  is 


Oi)'  *  K13  +  \ 


^131 

L  K133 

U 

-(ABe-<*K02)2 

I  ..  ..  /iT  ' 

L  Cos'1 1 

ABe-<*K02)V 

+  [K|3IK132  +  KI34  ^2  ' 

In 

H4 ' 

2  JJ 

*  K 


134 


K02  Cos'1  (  ABe‘?  ♦  K02  )  [cos'1  (  \Be"?  ♦  K02  )J 


V 1  -  K022 


[cos'^Be^  +  KoJ]4 


8 


Rewriting  the  inner  variable,  f; ,  in  terms  of  the  outer  variable,  h,  uses 
the  definition  ^  -  h/e .  Substituting  this  in  the  above  equation  gives 

K. 


(Ii),”Ki3  +  r  1 


"131 


133 


132  ^e"h/<*IC02 


+  K  02)1 


r  hx  \ 

L  Cos',(ABe'h/<.K()2)’| 

[K131K132  *  KI34  '  K15j 

In 

tan 

14 +  2  JJ 
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Rewriting  the  inner  variable,  f, ,  in  terms  of  the  outer  variable,  h,  uses 
the  definition  I;  -  h/e .  Substituting  this  in  the  above  equation  gives 


Since  eip(-l/e)  decreases  to  zero  rapidly  as  e-»0,  e"h/<  -  0  (Nayfeh, 
1981:260).  Therefore 

\  Be-h/<  ♦  K02  -  K02 

Thus,  the  outer  expansion  of  the  e1  inner  expansion  solution  for  a  is 

l(ai),]°“Ki5'Sin"1^02)  +  T^C2  ^  +  ,n  (494) 

V 1  "^02  L  M3  J  ■ 
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^02  ^  H4  ^OS  ^  ^+^02) 

Vl-^022 


[cos'1 

W<.k02) 

]3  3[cos'1(\Be^  +  K02)]5 

- 

6 

+  40  J 

[cos'1! 

W5*k02) 

]2  [cos'1(\Be‘^K02)]4 

2 

*  8  JJ 

Rewriting  the  inner  variable,  ^ ,  in  terms  of  the  outer  variable,  h,  uses 
the  definition  ^  -  h/e .  Substituting  this  in  the  above  equation  gives 
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Since  exp(-l/e)  decreases  to  zero  rapidly  as  e-*0,  e'h/€  -  0  (Nayfeh, 
1981:260).  Therefore 

X  Be  **  *  +  IC02  *  K02 

Thus,  the  outer  expansion  of  the  e1  inner  expansion  solution  for  Q  is 


now  used  to  correlate  the  inner/outer  expansions  and  outer/inner 
expansions.  From  Eq  (4.64)  the  matching  principle  states,  where  Y  is  an 
arbitrary  expansion 

[(v)°M(y)T 
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Thus,  the  inner/outer  expansions  and  outer/inner  expansions  derived 
above  are  equated,  resulting  in  the  original  inner  expansion  constants  of 
integration,  Kn,  being  expressed  in  terms  of  the  original  outer  expansion 
constants  of  integration,  Cn.  These  outer  expansion  constants  of  integration 
are  derived  from  the  initial  conditions  of  the  planetary  entry  problem  being 
examined. 


Blending  a.  Matching  Expansions.  Applying  the  matching 

principle,  Eq  (4.64),  to  the  inner/outer  and  outer/inner  expansions,  Eqs 
(4.87)  and  (4.92),  for  q,  gives 

*  C12 

3  ■  KI2IK02  ?  *  K!2 
(2  -  C0,C02)2 


Solving  for  Ki2  gives 


(4.96) 


Blending  U|  Matching  Expansions.  Applying  the  matching 

principle  to  the  inner/outer  and  outer/inner  expansions,  Eqs  (4.86)  and 
(4.91),  for  Uj  gives 


H 


,02  «P 


(  2 


^Cos'OCoj) 


1C, ,  +  2K 


1) 


Sin-’dC^)  ♦ 


2KftlIC 


AVI  -  K, 


2  L 


^12 ^  +  Kj2i  M 


Md-ic, 


JU¬ 


AB 


02 
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Solving  for  Kn  gives 


!iQLCll 
C, 


2K 


01 

2^01^02 


o,  j^*K12,  (C.  D 

VK02 


.2tC^Sin-i(Ko2) 


4(1 -Kq/) 


2\\  1 


yi- 


“02 


=[Kl2^Kmln-^ 


(4.97) 


Blending  a!  Matching  Expansions.  Applying  the  matching 

principle  to  the  inner/outer  and  outer/inner  expansions,  Eqs  (4.90)  and 
(4.94),  for  a,  gives 


CI2  (C01  +  0  +J^r(CQlC02  '  1 

T  7TT7  tT  *  cis  ‘  Ku '  sin'(K' 

V  1  "  C01  C02  J  v  2  *  C01(C02  +  1 ) 


Solving  for  K,j  gives 


Kn  -  C,«  + 


CI2  (C01  +  1)+^^C01C02^1) 


55  15  ( 1  -  C012C02)V2-C01(C02+1) 


+  Sin_1(K02) 


(4.98) 


Blending  I|  Matching  Expansions.  Applying  the  matching 

principle  to  the  inr-r/outer  and  outer/inner  expansions.  Eqs  (4.88)  and 
(4.93),  for  I,  gives 
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Solving  for  IC,3  gives 


(4.99) 


Applying  the  matching 


principle  to  the  inner/outer  and  outer/inner  expansions,  Eqs  (4.89)  and 
(4.95),  for  Q,  gives 


75 


As  in  the  derivation  of  e°  solutions  to  the  equations  of  motion,  the 
composite  expansion  or  solution  is  used  to  blend  the  distinct,  but 
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overlapping,  outer  and  inner  expansions  to  give  an  expansion  valid  over  the 
entire  altitude  domain.  The  composite  expansion  for  el  solutions  is  defined 
as  (Nayfeh,  1981:277) 


w, 


i  f  i  1°  o  i  r 

+  yj  -[Y|J  -  Yj  ^  Yj  -[YjJ 


(4.101) 


Also,  as  in  the  e°  derivation, 


is  used,  unless  otherwise  noted, 


since  it  typically  is  a  simpler  and  more  compact  expression  than 


u1  Composite  Expansion.  Applying  the  composite  expansion 

definition,  Eq  (4.80),  to  the  outer,  inner  and  inner/outer  expansions,  Eqs 
(4.27),  (4.59)  and  (4.86)  respectively,  for  Uj  gives 


“I  -  -C,  I]7h  *  ( *Be'h/<  *  K02  )2  exp[-  f  Cos’1  ( ♦  K„2  ) J 


where  K0I,  K02, 1C,2  and  K,21  are  given  by  Eqs  (4.76),  (4.75),  (4.96) 
and  (E.26),  respectively. 
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fli -Composite  Expansion.  Applying  the  composite  expansion  definition 
to  the  outer,  inner  and  outer/inner  (here  this  expansion  is  simpler  than  the 
inner/outer  expansion)  expansions,  Eqs  (4.28),  (4.60)  and  (4.92) 
respectively,  for  qt  gives 


j  11  +cl2'/n~h 

VC01(l  +  h)  12 

3 

(2-C01C02(l  +  h))2 


-ABe~h/c 


(4.103) 


Ij  Composite  Expansion.  Applying  the  composite  expansion  definition, 

to  the  outer,  inner  and  inner/outer  expansions.  Eqs  (4.29),  (4.61)  and  (4.88) 
respectively,  for  I,  gives 


•13 +  a 


K131KI32  ‘  KI3« 


(JL  v  V 

2  "  RI5 


A l 


lnl  tanl^  + 


n  Cos'1  (  ABe~h/e  +  K02 ) 


+  K 
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IC02  Cos'1  (\Be‘h/<  •  K02)  [cos'l(AB<fh/<  ♦  K02)]! 


VT 


"02 


Cos'^ABe-*1^  ♦  K02)]4 


8 


j  j 


(4.104) 


where  K02,  K ,5.  K,31,  KI32,  K,33,  and  ICI34  are  given  by  Eqs  (4.75), 
(4.98)  and  (E.30),  respectively. 
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a,  Composite  Expansion.  Applying  the  composite  expansion  definition 

to  the  outer,  inner  and  inner/outer  expansions,  Eqs  (4.31),  (4.63)  and  (4.90) 
respectively,  for  a,  gives 

(C|2  *  ChC02Vc<m  )(l+h)  -  (~r^+  C0ICI2 

C _ V  VL0» _ L  „ 

a|  *  /  ?  \  I  0  +  ^l5 

(l  “  C0|  C02 )  v‘^o, +  2(l+h)  -  C0IC02(l+hr 
-  Sin'1  (  NBe‘h/e ♦  K02 )  ♦  "T^-J 

V ) "^02 


where  K02  and  K15  are  given  by  Eqs  (4.75)  and  (4.98). 


0)  Composite  Expansion.  Applying  the  composite  expansion 

definition  to  the  outer,  inner  and  inner/outer  expansions,  Eqs  (4.30),  (4.62) 
and  (4.89)  respectively,  for  ft,  gives 
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♦  KH3M  tan 


L  Cos’^AB^^Koj)' 


14 


♦  K,aa  ^ 


4  ln(\Be~h/€*K02)  (NBe‘h/MC02)2" 


"in 

3 

k02k 

144  CoS 

1(\Be‘h/€>K 

02) 

+ 

V 1  -K022 

+  KI45 

[cos_1l 

Wh/'.K02) 

]3  3[cos'1 

(*B<>-h/‘.K02)]5 

6 

1  + 

40 

**146 

[cos'1! 

W^'.ICoj) 

J  [cos_1( 

\Be-h/<*K02)]4| 

2 

•f  ‘ 

8  JJ 

(4.106) 


where  JCq2,  Km,  ICj^2,  K^g,  K^,  K^,  and  K^g  are  given  by  Eqs 
(4.75)  and  (E32). 


Equations  of  Motion 


Once  the  composite  expansions  for  the  five  dependent  variables  (u,  q, 
I.  Q  and  a)  are  determined,  the  approximate,  analytical  solution  to  the  five 
equations  of  motion  are  given  as 

uc  -  u|j  ♦  (uj)e  ♦  CXe2) 

■  (-C<"TTh‘  K0,(^B«-h/‘*K02)2eip[-f  Cos-'^Be-h/'.K02)] 
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where  the  inner  expansion  constants  K01-K05  are  given  in  Eqs  (4.75M4.79) 
and  Kn-K15  are  given  in  Eqs  (4.96M4.100). 

The  above  set  of  five  equations  form  an  approximate,  analytical 
solution  to  the  set  of  five,  coupled,  first  order,  non-linear  ODEs  (Eqs  (4.1)- 
(4.5))  that  describe  the  trajectory  of  a  vehicle  entering  a  non-rotating 
planetary  atmosphere.  This  solution  set  is  first  order  accurate  with  respect 
to  the  small  perturbation  parameter  l/f)rs,  which  is  approximately  1/900 
for  earth  entry.  Since  a  solution  is  defined  as  Y  -  Y0  +  Y,e  +  0(e  )  and  the 
orders  of  Y0  and  Y,  are  simitar  by  definition,  addition  of  the  first  order 
solution  Y,e  increases  the  accuracy  of  the  solution  by  three  orders  of 
magnitude. 

Eqs  (4.107M4.1 1 1 )  are  relatively  simple  and  accurate  analytical 
solutions  to  the  complex  physical  system  of  atmospheric  entry.  They 
provide  a  readily  available  analysis  which  retains  the  subtleties  of  the 
original  system.  It  is  this  characteristic  of  analytical  solutions  which  makes 
them  more  favorable  than  obtaining  solutions  using  numerical  analyses. 
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V.  Validity  of  the  Solutions  to  the  Equations  of  Motion  and 
Comparison  to  Numerical  Solutions 


In  Section  IV,  analytical  solutions  are  derived  solving  the  set  of  five, 
coupled,  first  order,  non-linear  ODEs  that  describe  the  trajectory  of  a 
vehicle  entering  a  non-rotating  planetary  atmosphere.  These  solutions, 
given  as  Eqs  (4.107M4.1 1 1 ),  are  first  order  accurate  to  the  small 
perturbation  parameter  l/pr5,  which  is  approximately  1/900  for  earth 

entry.  They  provide  a  relatively  simple  and  accurate  solution  to  the 
complex,  non-linear  physical  phenomena  of  atmospheric  entry,  retaining  the 
trends  and  subtleties  that  are  lost  in  an  analysis  using  numerical  methods. 

Comparison  of  Analytical  and  Numerical  Results 

To  demonstrate  the  first  order  accurate  analytical  solution  derived  in 
Section  IV  (Eqs  (4.106M4.1 11)),  the  differential  equations  of  motion  (Eqs 
(4.1M4.5))  are  numerically  integrated  and  compared  to  the  results  from  the 
derived  analytical  solutions  for  the  same  set  of  initial  conditions. 

Numerical  Approach.  To  numerically  integrate  the  five  equations  of 
motions  derived  in  Section  III  (Eqs  (3.60M3.64)),  a  fourth-order  predictor- 
corrector  integrating  algorithm  is  used  (Wiesel,  1989:1 19-123).  The 
integrating  step,  Ah,  is  decreased  and  the  equations  of  motion  are 
repeatedly  integrated  until  the  resulting  data  from  the  above  integration 
remains  unchanged.  The  numerical  integrator  given  by  Wiesel  assumes  an 
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independent  variable  that  is  monotonically  increasing  or  decreasing.  For 
this  study,  the  non-dimensional  altitude  h  is  the  independent  variable,  but  it 
is  not  monotonically  decreasing,  since  the  skipping  of  the  re-entry  vehicle 
results  in  a  local  oscillation  in  altitude.  Thus,  the  numerical  integrator  is 
slightly  modified  to  account  for  this  altitude  oscillation  by  changing  the  sign 
of  the  altitude  increment  based  on  the  sign  of  the  flight  path  angle.(Ah  >  0  if 
y  >  0  and  Ah  <  0  if  y  <  0) 

Solution  Comparison.  To  compare  the  analytical  and  numerical 
techniques  to  solve  the  differential  equations  of  motion,  the  initial  conditions 
for  an  Apollo-type  reentry  vehicle  are  used  (Hillje,  1969:2-10).  Figures  FI - 
F5  (Appendix  F)  show  that  the  analytical  solution  give  very  accurate 
solutions  to  the  equations  of  motion.  The  advantage  of  having  a  first  order 
accurate  solution  is  evident  in  the  increased  solution  accuracy  around  skip 
points,  where  the  flight  path  angle  changes  sign,  and  at  low  altitudes,  where 
aerodynamic  forces  dominate. 

Derivation  Assumptions 

To  derive  analytical  solutions  to  the  five  ODE  equations  of  motion, 

Eqs  (4.1  M4.5),  several  assumptions  or  approximations  were  made  to  arrive 
at  analytical  solutions.  The  first  approximations  used  were  introduced 
when  the  small  perturbation  expansions  for  sin(a  +  be  ♦  ce  +0(e  ))  and 
cos(a  ♦  be  ♦  ce2  >0(e3))  were  performed.  Although  the  results  of  the 
expansions  were  not  exact,  they  were  correct  to  order  e  .  Since  the  final 
solutions  were  order  e°  and  e1 ,  this  assumption  does  not  effect  the  zero  and 
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first  order  solutions.  This  conclusion  is  also  applicable  to  the  binomial 
expansions  used  in  initial  expansion  of  the  equations  of  motion  and  the  zero 
and  first  order  matching  expansions. 

The  next  approximation  used  was  to  determine  the  small  perturbation 
expansion  for  tan(a  +  be  +  ce2  ♦0(e3)).  Although  the  expansion  is  correct  to 
order  e  ,  like  the  sine  and  cosine  expansions  above,  it  assumed  tan(a)  « 

1  /tan(be  +  ce2  ♦  0(e3)).  Since  b  and  c  »*  0(  1 )  and  e  <-  0(0.00 1 ),  this  is  a  valid 
assumption  for  most  common  orbits.  The  exceptions  to  this  inequality  are 
discussed  in  the  next  section. 

The  remaining  assumptions  used  during  the  derivation  of  the  zero 
and  first  order  solutions  were  made  to  solve  the  first  order,  inner 
expansions.  Through  numerical  analysis,  the  only  significant  approximation 
made  was  in  deriving  the  qt  inner  solution,  where 

♦  K02)eip[^Cos-'(^B<T?  ♦  K02)]  d? 

-  j(*Be'?»  K02)eipJ|cos"'(K02)]  d?  (5.1) 

Since  a  q,'  term  is  present  in  ail  of  the  first  order  ODEs,  the  error 
induced  by  this  approximation  is  present  in  all  of  the  first  order  solutions. 
This  is  reflected  in  Figures  E3-E6. 


Solution  Validity  and  Restrictions 

Even  before  the  above  assumptions  and  approximations  were  used, 
the  solutions  to  the  five  equations  of  motion  were  restricted  due  to 
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singularities  found  in  the  ODEs.  These  singularities  were  present  in  every 
ODE  except  dq/dh  and  involved  terms  of  sec(y),  csc(y)  and  csc(I).  Thus, 
the  solutions  to  the  ODEs,  whether  numerical  or  analytical,  become 
numerically  unstable  at  extremely  shallow  flight  path  angles,  lyl  «  1*  or  at 
extremely  steep  flight  path  angles,  y  <-  -90*.  The  latter  entry  trajectory  is 
not  realistic  since  it  results  in  enormous  aerodynamic  heating  and 
decelerating  forces.  The  solutions  also  encounter  instabilities  when  the 
entry  body  is  in  the  equatorial  plane,  I  -  0*. 

In  using  the  small  perturbation  expansions  for  the  tangent  of  I  and  y, 
the  approximations  impose  restrictions  on  the  validity  of  the  derived 
analytical  solutions.  As  mentioned  above,  the  expansions  assumed  tan(I)  « 
1/tandje  +  I2e2  +  0(e3))  and  tan(y)  «  I/tan(y,e  +  y2e2  +  0(e3)).  The 
assumptions  become  invalid  when  0(tan(I))  -  0(l/tan(I,e  +  I2e  +  0(e  )))  or 
0(tan(y))  -  0(  l/tan(y,e  ♦  y2e2  +  0(e3))).  For  earth  re-entry,  where  e  - 
1/900,  the  above  assumptions  break  down  when  85*  <  III  (or  lyl)  <  95’  For 
martian  re-entry,  where  the  mean  planetary  radius  is  smaller  than  the 
earth  s  and  the  atmosphere  is  thinner  than  earth's  (e  -  1/350  (Vinh  and 
others,  1980:5)),  the  above  assumptions  break  down  when  80’  <  III  (or  lyl)  s 
90*.  Thus,  the  approximate  solutions  become  invalid  at  extremely  steep 
entry  trajectories  or  near  polar  orbits. 

The  assumption  given  as  Eq  (5.1 )  was  used  to  facilitate  the  first 
order,  inner  solution  for  q(.  The  approximation  ABe"*  -  0  was  found 

accurate  until  the  entry  vehicle  penetrated  the  lowest  regions  of  the 
atmosphere,  when  the  approximation  slightly  underestimated  the  changes 
in  q,.  since  as  i^->0,  e'^  1.  Since  this  approximation  is  carried  forward  to 
the  first  order,  inner  expansion  ODEs  for  Ij,  Qj  and  a,,  in  the  form  of  the  qt 
terms  present  in  each  of  the  ODEs,  this  underestimation  trend  is  present  in 
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all  of  the  above  solutions,  as  shown  in  Figures  E4-E6.  The  induced  errors 
occur  at  very  low  altitudes,  where  terminal  course  corrections,  such  as 
lift/drag  modulation,  are  performed.  But  lift/drag  modulation  entails 
changing  the  lift/drag  coefficient,  which  was  assumed  constant  for  this 
study.  Thus,  the  errors  induced  by  the  above  approximation  are  de- 
emphasized. 
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VI.  Conclusions  and  Recommendations 


Conclusions 

Using  the  Method  of  Matched  Asymptotic  Expansions,  this  study 
developed  first  order  accurate,  analytical  solutions  (Eqs  (4.107M4.1 1 1 ))  to 
the  five,  coupled,  non-linear  equations  of  motion  describing  three- 
dimensional,  planetary  atmospheric  entry.  A  non-rotating  planet  and 
atmosphere  were  assumed,  as  well  as  a  constant  lift-to-drag  ratio  and 
ballistic  coefficient.  The  validity  of  the  developed  solutions  are  coordinate 
dependent  since  singularities  are  present  in  the  original  equations  of  motion. 
As  a  result  of  this  study,  the  following  conclusions  are  made 

1.  A  computerized,  symbolic,  algebraic  manipulator  greatly  reduced 
the  workload  in  generating  the  zero  and  first  order  asymptotic 
expansions.  Application  of  a  similar  program  could  be  used  in 
generating  similar  expansions  for  any  arbitrary  set  of  differential 
equations. 

2.  Exact  solutions  were  found  for  both  of  the  zero  order,  outer  and 
inner  asymptotic  expansions  and  the  first  order  outer  asymptotic 
expansions.  To  derive  analytical  solutions  for  the  first  order  inner 
asymptotic  expansions,  approximations  were  used  resulting  in 
underestimation  of  the  aerodynamic  turning  at  low  attitudes.  This 
inaccuracy  becomes  evident  at  low  altitudes  where  terminal 
maneuvers  are  initiated. 
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Recommendations 

Based  on  the  analysis  of  the  assumptions  and  the  results  of  this  study, 
the  following  recommendations  for  further  study  are  proposed. 

1.  Using  Mathematical  greatly  reduced  the  workload  of  finding  the 
outer  and  inner  asymptotic  expansions  to  the  equations  of  motion. 
Development  of  a  generalized  version  of  the  Mathematical 
program  used  in  this  study  would  result  in  a  program  generating 
the  nth  order  asymptotic  expansions  for  any  arbitrary  differential 
equation. 

2.  Due  to  approximations  in  small  perturbations  expansions,  the 
application  of  this  study  is  not  valid  for  entry  trajectories  near 
polar  orbits.  A  different  set  of  coordinate  transformations  and 
subsequent  expansions  should  be  applied  to  extend  the  domain  of 
solutions  near  this  orbital  state. 

3.  Due  to  approximations  in  the  first  order,  inner  expansion  solution 
for  q,  the  flight  path  angle,  the  solutions  found  in  this  study 
become  inaccurate  for  low  vehicle  altitudes.  Further  investigation 
should  be  undertaken  to  find  either  exact  solutions  to  the  first 
order,  inner  expansion  ODEs  or  better  approximations  which 
would  increase  the  accuracy  of  the  analytical  solutions  at  low 
altitudes. 
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Appendix  A: 

Transformation  of  Planetocentric  Coordinates  to  Orbital 
Elements  Using  Spherical  Trigonometry 


This  appendix  derives  the  relations  between  the  planetocentric 
angles  from  the  equations  of  motion(6,  4>  and  (j>)  and  the  classical  orbital 
elements  (a,  Q  and  1).  Since  this  study  models  the  planet  as  a  sphere  (Vinh 
and  others,  1980:2),  spherical  trigonometry  is  used  to  derive  relationships 
between  the  planetocentric  angles  and  the  orbital  elements  (Bain.  1989). 

Figure  A1  shows  the  relationship  between  the  two  related  variable 
sets  in  an  osculating  orbit  (Vinh  and  others,  1980:257).  For  clarity.  Figure 
A2  shows  the  pertinent  variables  as  a  spherical  triangle  formed  by  the  arcs 
of  three  great  circles  on  the  surface  of  a  sphere  (planet).  By  definition,  the 
interior  angles  are  angles  between  the  curved  line  segments  and  the 
exterior  angles  are  the  angles  between  the  linear  segments  emanating  from 
the  origin,  point  0.  The  following  is  a  list  of  the  interior  and  exterior  angles 
used  in  the  subsequent  derivation. 

Eiterior  Angles:  a-e-Q  b-<j>  c  -  a 

Interior  Angles:  A  -  tt/2  -  tj>  B  -  I  C  -  tt/2 

To  relate  the  above  two  variable  sets,  the  two  fundamental  formulas 
of  spherical  trigonometry,  the  law  of  sines  and  the  law  of  cosines,  are  used 
(Fitzpatrick,  1970:1 18-120).  While  the  law  of  sines  relates  the  ratios  of 
interior  and  opposing  exterior  angles,  the  law  of  cosines  relates  the  cosine  of 
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Figure  A 1 .  Reference  Coordinate  System  and  Orbital 

Elements 


an  interior  angle  as  the  sum  of  the  products  of  the  cosines  of  the  other  two 
interior  angles  and  the  products  of  the  sines  of  the  other  two  interior  angles 
with  the  cosine  of  the  opposing  exterior  angle. 


Law  of  Sines: 
Law  of  Cosines: 


sin(a)  sin(b)  sin(c) 
sin(A)  "  sin(B) "  sin(C) 

cos(a)  -  cos(b)  cos(c)  ♦  sin(b)  sin(c)  cos(A) 


(A.l) 

(A.2) 


To  derive  the  three  spherical  trigonometric  relationships,  known 


values  are  first  substituted  into  the  law  of  sines 


sin(6  -  fl)  sin(4>)  sin(a) 
sin(n/2  -  <|>)  "  sin(I)  "  sin(n/2) 


(A.3) 
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Figure  A2.  Spherical  Triangle  Relating 

Planetocentric  Angles  to  Orbital  Elements 


Since  sin(rr/2  -  4> )  -  cos(<|> )  and  sin(n/2)  -  1,  two  relations  arise  from 
Eq  (A.3)  and  making  use  of  the  common  trigonometric  relationship, 
sin2(x)  ♦  cos2(x)  -  1: 


sin(a) 


sin(4>) 

sin(I) 


or 


2 

cos  (a)  -  1  - 


sin2(<|>) 

sin2(I) 


(A.4) 


2/«  ^  .  sinz(d>)  cosz(0> ) 

or  cos  (0  -  fi)  -  1  - - T  ■*  -  -Y-~ 

sin‘(I) 


(A.5) 


Using  the  law  of  cosines  gives  the  first  relationship  as 

cos(a)-  cos(6  -  0)  cos(<J>)  +  sin(0  -  fi)  sin(<f>)  cos(tt/2)  or 

cos(a)-  cos(0  -  0)  cos(4>)  (A.6) 
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Squaring  (A.6)  and  substituting  Eqs  (A.4)  and  (A.5)  into  (A.6)  gives 


sin2(<|>)  (  sin2(<t>)  cos' 

sin2(I)  sin2(I) 


cos2(4>) 


Simplifying  the  above  expression  gives  the  second  relationship: 


cos(I)  -  ±cos(<J>)  cos(4>) 


(A.7) 


To  derive  the  third  relationship,  Eq  (A.3)  also  gives 


sin(6  -  ft)  - 


sin(I) 


(A.8) 


But  from  Eq  (A.7) 

cos(I) 

cos,,|,)-^ 


(A.9) 


Substituting  Eq  (A.9)  into  Eq  (A.8)  gives  the  third  relationship: 


sin(9  -  ft) 


tan(<j>) 

tan(I) 


(A.10) 


Thus  Eqs  (A.6),  (A.7)  and  (A.  10)  give  the  three  spherical 
trigonometric  relationships  required  to  transform  the  pianetocentric  angles 
to  classical  orbital  elements.  For  these  equations  to  be  useful,  they  will  be 
reformed  to  match  the  equations  of  motion  Eqs  (3.46M3.50).  Eq  (A.4)  is 
already  in  the  correct  format.  Repeating  Eq  (A.6) 


cos(e  -  ft)  - 


cos(q) 

cos(4>) 


(All) 


But  from  the  law  of  cosines 
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(A. 12) 


cos(0  -  Q)  -  cos(<f>)  cos(a)  +  sin(<J>)  sin(a)  cos(tt/2  -  <|> ) 
-  cos(4>)  cos(a)  +  sin(<J>)  sin(a)  sin(<j>) 


Equating  Eqs  (A.l  1)  and  (A.  12)  and  simplifying  gives 
cos(a)  sin(I) 


sin(tM 


cos(4») 


Substituting  Eq  (A.4)  into  Eq  (A.  13)  gives 
.  , .  v  tan(4>) 


(A.13) 


(A. 14) 


Again  taking  the  law  of  cosines 

cos(rr/2  -  <J>  )«  cos(I)  cos(tt/2)  +  sin(I)  sin(rr/2)  cos(9  -  Q)  or 

sin(4>)  -  sin(I)  cos(0  -  Q)  (A.15) 

To  summarize,  the  three  transformation  relationships  between 
planetocentric  angles  (4>,  <J>  and  6)  and  classical  orbital  elements  (a,  Q  and  I) 
are: 


sin(<|>)  -  sin(a)  sin(I) 

(A.4) 

cos(I)  -  cos(<l>)  cos(4>) 

(A.7) 

sin(4>)  -  sin(I)  cos(0  -  Q) 

(A.15) 
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Appendix  B: 

Mathematical  Code  IQ  Generate  Zero.  First,  and  Second  Order 

Asymptotic  Expansions 


This  appendix  lists  the  Mathematics™  code  used  to  generate  the  zero, 
first  and  second  order  outer  and  inner  expansions  to  the  five  equations  of 
motion 


Program  Structure 

The  Mathematica™  program  is  structured  to  input  a  matrix  whose 
first  column  is  composed  of  the  number  of  ODEs,  followed  by  the  ODEs 
themselves.  The  program  substitutes  the  small  perturbation  expansions  for 
the  dependent  variables,  Eq  (4.6)  and  also  for  trigonometric  and  algebraic 
functions  of  the  dependent  variables.  For  the  outer  expansions,  the  program 
multiplies  out  all  of  the  terms  composed  of  sums  of  constant  coefficients  and 
powers  of  c.  The  program  then  combines  the  coefficients  of  e°,  e1  and  e* . 
The  6°  terms  are  the  zero  order  outer  expansion  ODEs,  el  terms  are  the  first 
order  outer  expansion  ODEs  and  e  terms  are  the  second  order  outer 
expansion  ODEs.  To  perform  the  inner  expansions,  the  independent 
variable  h  is  changed  to  the  magnified  variable  ^  by  the  stretching 
transformation  h  -  .  The  program  then  multiplies  and  collects  terms  as 

in  the  outer  expansions  just  completed.  The  e°  terms  are  the  zero  order 
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inner  expansion  ODEs,  e1  terms  are  the  first  order  inner  expansion  ODEs 
and  e2  terms  are  the  second  order  inner  expansion  ODEs. 

The  core  of  the  program  is  the  sequence  which  multiplies  out  the 
small  perturbation  expansions  of  the  dependent  and  independent  variables 
and  collects  coefficients  of  powers  of  e.  These  steps  are  executed  by  the 
Mathematica™  functions  Expandl)  and  Coefficient!).  Expand!)  writes  the 
products  of  polynomials,  in  powers  of  c,  as  a  simple  sum  of  terms  of 
constants  coefficients  and  powers  of  e  with  all  products  expanded  out. 
Coefficient!)  collects  coefficients  of  a  prescribed  power  of  e  from  the  the 
above  sum  of  terms  (Wolfram,  1988:381-384). 


Program  Listing 

*X*S**(*SXt**t*t*****t*X**>******t**S*S(«S**tt**(tS«*t*S*****t*SS******XI 

Expansions : :  usage  -  "This  Mathematica  function  gives  the  zero,  first 
and  second  outer  and  inner  asymptotic  expansions  for  a  set  of  first 
order  differential  equations  of  motion  (EOM).  The  EOM  are 
inputted  via  a  matrix  called  ode*,  ode!!  1,1)1  is  the  number  of  EOM 
being  expanded.  ode!l2,l)l,  ode({3,I)l,  etc  are  the  actual  EOM.  The 
expansion  parameter  must  be  called  e'  for  the  program  to  work." 


ExpansionslodeJ :  - 

(*  Defining  small  perturbation  expansions  needed  for  outer  and  inner 
expansions  *) 


(*  Defining  the  small  perturbation  expansions  for  the  dependent 
variables  in  terms  of  the  small  parameter,  e  *) 
u  -  uO  +  ul*e  ♦  u2*e‘2  ; 

q  -  qO  *ql*e  ♦  q2*e‘2  ; 
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i 

*  iO 

+  iPe 

+  i2*e‘2  ; 

0 

-  oO 

+  ol*e 

+  o2*e‘2  ; 

a 

-aO 

+  aPe 

♦  a2*e*2 ; 

g 

-  gO 

♦  gPe 

♦  g2*e*2  ; 

(*  Defining  the  small  perturbation  expansion  for  trigonometric 
functions  *) 

Sinlx0_  +  xl_*e_  +  x2_*e_‘2)  -  Sin  1x0)  +  xPCos[xO)*e  +  ((x2*Cos[x01)  - 
(xr2*sin(x0)/2))*e‘2 ; 

Coslx0_  ♦  xl_*e_  ♦  x2_*e_‘21  -  CoslxO]  -  xl*Sinlx01*e  -  ((x2*Sin[x01)  + 
(xP2*cos(xO)/2))V2  ; 

Tan|x0_  +  xl_*e_  +  x2_*e_‘2)  -  TanlxO]  +  (xl/CoslxO)'2)*e  + 
((x2/Cos[x0r2)  ♦  (x  1  ‘2*sin(x0 )/cos(x0  )‘3  ))*e‘2  ; 


(*  Defining  the  small  perturbation  expansion  for  the  reciprocal  of 
trigonometric  functions  *) 

Sin  [x0_  ♦  xl_*e_  *  x2_*e_‘2M  -  (1/SinlxO))  - 

(x‘*CosIxO]/Sin|xOr2)*e  *  ((xl*2*(l  ♦  Cos[x0r2)/2*Sinlx0)*3)  - 
(x2*Cos|xO)/Sin|xO]‘2))V2 ; 

Coslx0_  +  xl_*e_  ♦  x2_*e_*2r-l  -  ( 1  /CoslxO))  + 

(x  1  *Sin|i01/Cos|i0]‘2 )*e  ♦  ((x  1*2*0  ♦  Sinlx0r2)/2*Coslx0)*3)  ♦ 
(x2*SinIx0)/Coslx0r2))*e‘2  ; 

Tan  [x0_  ♦  xl_*e_  ♦  i2_*e_7)*-l  -  (1/TanlxOI)  -  (xl/Sin|x0r2)*e  ♦ 
((ir2*Cos|x0r2/Sin[x0r3)  -  (x2/Sin|xO)‘2))*e‘2 ; 

(*  Defining  the  small  perturbation  expansion  for  algebraic  functions  *) 

(x0_  +  x l_*e_ ♦  x2_*e_‘2)‘- 1  -  (1/xO)  -  (xl/x0*2)*e  ♦  ((xl‘2/x0*3)  - 
(x2/i0‘2))*e‘2 ; 


(*  Defining  the  exponential  of  a  large  negative  number  is  zero,  or 
allowing  the  expansion  to  neglect  exponentially  small  terms  *) 

(E*(h/e)M  -  0  ; 
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(*  Performing  the  outer  expansions  of  the  EOM  *) 


(*  Multiplying  out  all  the  terms  in  the  EOM  *) 

Dolodelli ,  21)  -  Expandlodelli ,  1  ]]] ,  0 , 2 ,  ( 1  ♦  ode[[  1 , 1  ])))] ; 

i  - . ; 


(*  Grouping  the  coefficients  of  e  raised  to  the  0,1  and  2  powers.  The 
zero  order  outer  expansions  are  located  in  ode!12 , 311-ode(16 ,311, 
the  first  order  outer  expansions  are  located  in  ode[[2 ,4ll-odell6 , 
4)1,  and  the  second  order  outer  expansions  are  located  in  odej[2 , 
5H-ode[[6.5H.  *1 

Dolodelli, 311  -  Coefficientlodelli,  211, e,  01,  {i,  2,(1  +  odelll .  Ill))) 

i  - . ; 

Dolodelli, 411  -  Coefficientlodelli , 211, e,  1] . (i . 2,(1  +  odelll  ,111)11 

i  - . ; 

Dolodelli ,  511  -  Coefficientlodelli , 211 . e ,  21 , (i , 2 , ( 1  +  odell  1 . 1  ]])}] 

i  - . ; 


(*  Performing  the  inner  expansions  of  the  EOM,  which  transforms  h 
to  z  by  the  definition  h  -  z*e) 


(*  Multiplying  out  all  the  terms  in  the  EOM  *) 

Dolodelli ,  211  -  e’odelli .  1 1) ,  {i ,  2 ,  ( 1  ♦  odell  1 , 1  U))l ; 

i  -  • ; 

(*  Implementing  the  definition  relating  h  and  the  stretching  variable  z 
*1 

h  -  z*e ; 

(*  Multiplying  out  all  the  terms  in  the  EOM  *) 

Dolodelli ,  211  -  Expandlodelli ,  2111 ,  (i ,  2 ,  ( 1  ♦  odell  1 , 1 11)11 ; 
i  -  . ; 
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(*  Grouping  the  coefficients  of  e  raised  to  the  0,1  and  2  powers.  The 
zero  order  inner  expansions  are  located  in  ode([2 , 61J-ode([6 , 611, 
the  first  order  inner  expansions  are  located  in  ode([2 , 61]-ode[[6 , 
6]],  and  the  second  order  inner  expansions  are  located  in  ode{(2 , 
6ll-ode([6 , 61J.  *) 

Dolodelli ,  611  -  Coefficientlodedi , 211 . e , 01 ,  (i ,  2 ,  ( 1  ♦  odell  1 , 1  ]]))] 
i  - . ; 

Dolodelli ,  711  =  Coefficientlodelli ,  211 ,  e ,  1 1 ,  (i .  2 ,  ( 1  +  odell  1 , 1 11))1 
i  *  ; 

Dolodelli, 811  -  Coefficientlodelli , 211  ,e, 2) . (i , 2,(1  +  odell  1 , 111)11 
i  - . ; 

ytxtttmmmtmyttttmitmttttmmttttututmtttttttttmmstt 


Sample  Incut 


The  following  listing  is  a  sample  input  required  to  execute  the  above 
program.  It  is  same  the  five  ODEs  derived  in  Section  IV  and  describe  the 
flight  trajectory  of  a  lifting  body  entering  a  planetary  atmosphere. 


odell  1 , 1 11  -  5 

ode([2,  111  -  u/(l  +  h)  -  (2*b*u*(l  *  rTanlgl))/(E*(h/e)*e*Sin[gJ) 

odell3,  111  -  -((b*l)/(E‘(h/e)*e)-(q*(  1  -  q*2/u))/(l  *  h) 

ode((4, 11]  -  1  /(( 1  ♦  h)*Tanlgl)  -  (b*d*Sinlal)/(E‘(h/e)*e*Coslgl*Sin[gJ*Tanlil) 

odell5,  111  -  (b*d*Sinlal)/(E‘(h/e)*e*CoslgrSinlgl*Sinlil) 

odell6,l  ]]  -  (b*d*Coslal)/(E‘(h/e)*e*Coslgl*Sinlgl) 
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Appendix  C: 

Derivation  of  Frequently  Used  Tavlor  Series  and  Binomial 

Expansions 


This  appendix  derives  many  Taylor  series  and  binomial  expansions 
frequently  used  in  this  study. 


Small  Perturbation  Expansions  for  Trigonometric  Functions 

To  expand  the  sine,  cosine  and  tangent  functions  of  small 
perturbations  (powers  of  e)  into  linear  combinations  of  powers  of  the 
perturbations,  the  smail  angle  formulas  for  the  above  functions  are  needed. 
The  small  perturbation  here  is  e,  where  e  <  1 
Since  cos  (x)  is  defined  as 


cosu)  - 1  fcgat 

n-o  (2n)f 

,  .  ,  (ae)2  (ae)4 

cos(ae)  -  1  -  ~2\~  * "  41""" + 

-  1  ♦  |[  e2+  0(e4) 


Since  sin(x)  is  defined  as 


00 


2o*  I 


Sin(J)-z<- i£a 

0To  (2n+l)l 


(C.l) 
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To  derive  the  small  perturbation  expansions,  angle-sum  relations  are 
also  needed. 

cos(a  ♦  b)  -  cos(a)  cos(b)  -  sin(a)  sin(b)  (C.4) 


sin(a  +  b)  -  sin(a)  cos(b)  ♦  cos(a)  sin(b) 


tan(a  ♦  b) 


tan(a)  tan(b) 

1  -  tan(a)  tan(b) 


(C.5) 

(C.6) 


2  3 

In  deriving  the  small  perturbation  expansions  for  sin(a  +  be  ♦  ce  +de  ) 
and  cos(a  ♦  be  ♦  ce2  ♦  de3),  Eqs  (C.4)  and  (C.5)  are  repeatedly  used.  The 
following  example  illustrates  this  by  deriving  the  expansion  for 
sin(a  ♦  be  ♦  ce2  ♦  de3). 

sin(a  ♦  be  ♦  ce2  ♦  de3) 

2  3  2  3 

-  sin(a)  cos(be  ♦  ce  +  de  )  ♦  cos(a)  sin(be  +  ce  ♦  de  ) 
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-  sin(a)  (cos(be)cos(ce2  +  de3)  -  sin(be)sin(ce2  +  de3)) 

+  cos(a)  Isin(be)  cos(ce2  ♦  de3)  +  cos(be)  sin(ce2  +  dc3)! 

*  sin(a)  {cos  (be)lcos(ce2)  cos(de3)  -  sin(ce2)  sin(de3)l 
-  sin(be)  (sin(ce2)  cos(de3)  +  cos(ce2)  sin(de3)]} 

+  cos(a)  {sin(be)  (cos(ce2)  cos(de3)  -  sin(ce2)  sin(de3)l 
+  cos  (be)lsin(ce2)  cos(de3)  ♦  cos(ce2)  sin(de3)]} 


Substituting  in  Eqs  (C.l)  and  (C.2) 

2  o 

sin(a  +  be  +  ce  +  de  )  -  sin(a)  +  bcos(a)e 

+  |ccos(a)  -  —  s^n-a^e2  +  o(e3) 

Similarly.  cos(a  +  be  +  ce2  ♦  de3)  becomes 

cos(a  ♦  be  ♦  ce2  +  de3)  -  cos(a)  -  bsin(a)e 

.  ,  .  b2cos(a))  ?  qv 
-  csm(a)  ♦  — 2 — le  *  0(r) 


(C.7) 


(C.8) 


To  derive  the  expressions  for  the  reciprocals  of  the  above 
relationships,  binomial  coefficients  and  Pochhammer  symbols  are  used 
(Andrews,  1985:10-11,  273).  By  definition,  where  (-n)0  -  1  and 
(-n)k  -  (-n)  (-n  +l)...(-n  +  k  -  1) 


oo 


CO 


(a*  bf  - 

k«0 


bk 


(C.9) 


As  an  example  to  find  the  reciprocal  of  a  binomial  series 
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CO 

- 1 - -  y  M|a'l'k (be  +  ce2  +  de3f 

a  +  be  +  ce2  +  de3  k30  ^ 


-  Y  tiHlka-(l.lt)(be.C62.d€3f 

^  k! 

k*0 


a 


V 


e  ♦ 


a‘ 


a* 


e2  ♦  0(e3) 


(C.10) 


Substituting  Eqs  (C.7)  and  (C.8)  into  (C.10)  gives 


(sin(a  +  be  ♦  ce2  ♦  de3))"1 

1  bcos(a) 

fb2(  1  ♦  cos2(a))  ccos(a) 

'  sin(a) "  sin2(a)  * " 

(  2sin3(a) 

sin2(a) 

2  3  1 

(cos(a  +  be  +  ce  +  de  )) 

1  bsin(a) 

fb2(l  ♦  sin2(a))  csin(a)) 

’  cos(a) T  cos2(a)  *  ” 

[  2cos3(a) 

cos2(a)J 

e2  +  0(e3)  (C.ll) 


e2  +  0(e3)  (C.12) 


2  3 

To  derive  the  small  perturbation  expansion  for  tan(a  ♦  be  +  ce  ♦  de  ), 
Eq  (C.6)  can  be  used  repeatedly,  with  the  assumption  tan(a)  tan(be  ♦  ce2  ♦ 
de3)  -  0  (alternatively,  tan(a) «  l/tan(be  +  ce2  +  de3)).  This  assumption 
worsens  as  a-»n/2  radians  or  as  e  increases.  Due  to  this  approximation, 
Karasopoulos  limited  his  study  to  I  *  75°  (Karasopoulos,  1988:6)  To  avoid 
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this  approximation.  Eq  (C.7)  is  divided  by  Eq  (C.12),  since  tan(x)  - 
sin(x)/cos(x). 


?  .  sin(a  ♦  be  +  ce2  ♦  de3) 

tan(a  +  be  ♦  ce  +  de  )  -  — - - 5 - j- 

cosia  ♦  be  ♦  ce  +  de  ) 


sin(a)  +  bcos(a)e+  |ccos(a)  -  ~,s^a^  e2 


1  bsin(a)  fb2(  1  »  sin2(a))  csin(a) 


cos(a)  cos2(a) €  {  2cos3(a)  cos2(a)J 


♦  0(e3) 


Simplifying  and  grouping  terms  of  powers  of  e  gives 
tan(a  +  be  ♦  ce2  +  de3) 


/b2sin(a) 


tan(a)  +  j  e  +  i  +  a 

cos  (a)  cos  (a)  cos  (a) 


e2  +  0(e3)  (C.13) 


To  derive  the  reciprocal  of  the  above  expression.  Eq  (C.8)  is  divided 
by  Eq  (C.l  1).  since  l/tan(x)  -  cos(x)/sin(x). 

/.  /  .  2  .  3w-i  cos(a  +  be  +  ce2  +  de3) 

(tan(a  ♦  be  ♦  cc  *  de0))  1  -  — ; - 5 - 

sin(a  +  be  +  ce  ♦  de  ) 


-  Jcos(a)  -  bcos(a)e-  |csin(a)  +  ~ c^a^|e2 


1  bcos(a) 

b2(  1  ♦  cos2(a)) 

ccos(a) 

2 

sin(a)  sin2(a)  € 

k  2sin3(a) 

sin2(a); 

e 

Simplifying  and  grouping  terms  of  powers  of  e  gives 
(tan(a  ♦  be  ♦  ce2  +  de3))'1 
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I 


This  section  derives  many  of  the  frequently  used  binomial  expansions 
used  in  the  derivations  in  this  study.  Again  binomial  coefficients  and 
Pochhammer  symbols  are  used  and  are  defined  in  Eq  (C.9). 


oo 


va 


— V — *  Y  (?)  (be  *  c«2  * 

+  be  +  ce2  +  de3J  k„0  k 


co 


-  V  (-lyjlflfc  a-(o»k)  (be  +  ce2  +  de3 f 

k! 


k«0 

oo 


y  ♦  ce2  ♦  de3f 

tTo  ^ 


(C.15) 


For  example,  as  shown  above  in  Eq  (C.10),  for  n  -  1 
_ 1 _  (1)0  (be  ♦  ce2  +  d€3)0 


(a  *  be  ♦  cc2  +  dc3)1 


(0!)  a 


(l).(be  +  ce2  +  dc3)1  (l)2(be  ♦  ce4  +  de3)*  3 

(l!)a2  (21)  a3 

1  be  +  ce2  +  de3  2(be  *  ce2  ♦  de3)2  3, 

"  - 2 - ♦ - 71 - +  0<e  ) 


.2  . 


a 


a" 


2a 
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(C.16) 


€2  ♦  0(€3) 


As  another  example,  as  shown  above  in  Eq  (C.16).  for  n  -  Vi 


(It)  a2  (2!)  a2 


1  be  +  ce2  +  de3  3(be  +  ce2  +  de3)2  _ , 

"7  " - 3 - + - 5 - +  ) 

a2  2a2  8a2 

—  e2  +  0(e3)  (C.17) 

2a2 


-L  JL  (3b2 

1  “  3 €  +  5 

a2  2a2  t8a2 
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Appendix  D: 

Derivation  of  Solutions  for  Outer  Expansions  of  the  Equations 

of  Motion 


Methods  oLSoimton 

The  solutions  to  the  two  sets  of  five,  coupled,  first  order,  linear 
nonhomogeneous  ordinary  differential  equations  (ODEs)  found  in  this 
appendix  are  derived  below.  Three  methods  are  used  to  solve  these 
equations.  The  first  method  used  is  the  Method  of  Separation  of  Variables, 
where  the  ODE  and  its  solution  are  given  as  (Beyer,  1984:315) 

fr‘«y)  -»  J «y) dy * J f(l) dl - c  (D.l) 

C  is  the  constant  of  integration. 

The  second  method  used  to  solve  ODEs  found  in 
Method  of  Integrating  Factors,  where  the  ODE  has  the 

^+P(x)y-Q(x) 

The  integrating  factor  has  the  form 

u  -  e  Jp(,)  «•»  (D.3) 

and  the  solution  is  given  as  (Beyer,  1984:315 ) 
ye  Jp<,)  *  -  foil)  e  Ipl,)  *  ♦  C  or 


this  study  is  the 
form 

(D.2) 
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(D.4) 


v 


,!Q(xle  JP(xldx 

e  lp<») * 


*C 


The  third  and  final  method  is  a  variation  of 
known  as  Bernoulli's  Equation  (Beyer,  1984:315). 

■j*r «•  P(x)y  -  Q(x)yn 


the  second  method  and  is 
The  ODE  is  given  as 


(D.5) 


Again, the  integrating  factor  has  the  form 
.  e  Jp(x)  dx 


v 


(D.6) 


and  the  solution  is  given,  for  n  -  1,  as  (Beyer,  1984:315) 

y(l-n)e  (l-njptx)  dx  .  (1  .  n)/o(x)  e  ^  ")Jp(x)  dx  +  c 


or 


y  * 


( (1  -n)fo(i)ea'n)Jp(,)d,*cV~n 

e  (l-n^Jp(x)  dx 


V 


(D.7) 


Although  some  of  the  ODEs  are  coupled,  the  coupling  is  avoided  by 
solving  by  ODEs  in  a  judicious  order.  The  du/dh  ODE  is  solved  first, 
followed  by  the  dq/dh  ODE  is  second  and  finally  the  da/dh  ODE  is  solved. 
The  dl/dh  and  dQ/dh  ODEs  are  independent  of  the  other  ODEs  and  are 
solved  in  the  order  presented  for  consistency. 


Constants  of  Integration  Notation 
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As  shown  in  Section  IV,  there  are  ten  ODEs  expanded  from  the 
equations  of  motion.  Accompanying  the  solution  of  each  of  these  ODEs  is  a 
constant  of  integration.  To  simplify  the  bookkeeping  involved  in  defining 
all  these  constants,  the  following  notation  is  used  for  the  outer  expansion 
solutions: 

The  constant  of  integration  is  given  as  C^,  where  i  denotes  the  order 
(e1  --  e°  or  e1)  of  the  solution  and  j  denotes  the  variable  associated  with  this 
constant. 

j _ variable 

1  u 

2  q 

3  I 

4  Q 

5  a 

Thus  C0I  is  the  constant  of  integration  for  the  variable  u0  and  CH  is 
the  constant  integration  for  the  variable  Q,. 

Quter  Expansion  Terms 

The  differential  equations  corresponding  to  the  e°  equation  of  motion 
terms  were  derived  in  Section  IV  and  are  repeated  below: 


duQ 
dh  * 

_Ufi_ 

'  1  ♦  h 

(D.8) 

dq0 
dh  “ 

l*h(u0 

(D.9) 

dlfi 
dh  " 

0 

(D.10) 

d^ 
dh  “ 

0 

(Dll) 
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«2<l  m  . -  i 

dh  ( 1  +  h)  tan(y0) 


(D.12) 


du0/dh  Equation  Solution.  The  Method  of  Separation  of  Variables 
solves  Eq  (D.8).  Rearranging  Eq  (D.8)  into  form  of  Eq  (D.l)  gives 


du0  u0  au0  dh  fauft  f  dh 

dh  -  *i+h  “ST-TTh  Urdu+  TThdh-coi 


du0  dh 


Solving  the  above  ODE  gives 


ln(u0)  +  ln(  1  +  h)  -  C0I  or  u0-yj 


Solution:  u0  * 


(D.13) 


dq0/dh  Equation  Solution.  This  ODE  is  solved  by  substituting  the 
solution  of  u0,  Eq  (D.13),  into  Eq  (D.9).  Thus,  the  ODE  now  becomes 


dh  1  +  h  u 


1  ^  C0I 

\1  +  h 


Cq i  1  +  h 


The  ODE  now  has  the  form  of  Bernoulli’s  equation  Eq  (D.5),  where 


P(xl- 


1  +  h 


Q(x)  -  7T—  n  -  3 


From  Eq  (D.6) 
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u  -  e 


Jp(*)dx_e  .  in(i  .  h) 

Thus,  Eq  (D.7)  gives 

-2e  -21n(l  ♦  hj  ,  (,2)  f_L^  -21n(l  ♦  h 


Qo 


'01 


)dh  -  C02  or 


-2 

f  1 

C0I> 

(1  +h)2 

'02 


Integrating  the  right  hand  side  of  the  above  equation  gives 

1  2  -c 


(q0(  1  *  h))2  C0,(  1  ♦  h)  02 


Solving  for  q0  gives  the  desired  solution. 

I 


Solution:  q0 


1  +  h)  "  C02^1  +  h)2 
C01  V^02 


(D.14) 


dl0/dh  Equation  Solution.  The  solution  to  this  ODE  (Eq  (D.10))  is 

trivial. 

Solution:  I0-C03  (D.15) 

dfy/dh  Equation  Solution.  The  solution  to  this  ODE  (Eq  (Dll))  is 

trivial. 

Solution:  Qq-Cq^  (D.16) 
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daft/dh  Equation  Solution.  This  ODE  is  solved  by  first  rewriting  the 
term  tan(y0)  in  Eq  (D.12)  in  terms  of  q0,  where,  by  definition,  q0  -  cos(y0). 

1  1  1  1  1 

/X7T 

cos(^)  cosier  — ,  Vqo2 


But  from  Eq  (D.14) 

~K  *  'r~^  1  ♦  h)  -  C02(  1  ♦  h)2 

q0  Coi  02 


Substituting  the  above  two  relations  into  the  ODE,  Eq  (D.12)  gives 
dap  _ l  1 


dh"(l+h)  f 2  \ 

1 +  "  C02^  +  hr  - 1 


The  Method  of  Separation  of  Variables  is  used  to  solve  this  ODE. 
From  Eq  (D.l ),  the  solution  has  the  form 

1  1 


Jda0- 


(i  *  h)  n  \ 

1  *  h)  -  c02(i  ♦  hr  - 1 


dh  +  C, 


05 


The  above  integral  has  the  following  form  and  solution  (Beyer, 
1984:257), .where  x  -  1  ♦  h,  dx  -  dh,  a  -  -1,  b  -  2/C0l  and  c  -  -C02. 


r 


dx 


.\(  bx  ♦  2a  \ 
Ixl  Vb2-4acj 


♦  C 


From  the  definition  of  h  (Eq  (3.22)),  h  >  0,  therefore  |1  +  hi  -  1  +  h. 
Using  this  relation  and  substituting  the  above  relation  into  the  ODE  gives 
the  solution 
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Oq  *  r  Sin  ^ 


VI 


2(1  +h) 

Cp| 


-2 


(1  +  h) 


(-LS 2 


4C 


'01, 


02 


+  c 


05 


Simplifying  the  above  equation  gives 

C, 


Solution:  Oq  -  Sin 


-l 


1  - 


QL 


1  +  h 


LV 


1  ■  C0I  C02-| 


+  C, 


05 


(D.17) 


Outer  Expansion  e*  Terms 


The  differential  equations  corresponding  to  the  e1  equation  of  motion 
terms  were  derived  in  Section  IV  and  are  repeated  below: 


dqt 

dh 


dl 

dh 

dQ, 

dh 


u, 


1  ♦  h 


do 

1  ♦  h 


(2Mi  SaJkl 


Uf 


u 


0  ) 


_aj_ 

1  *  h 


(„  2  \ 

^-1 

Un 


0 

0 


(D.18) 

(D.19) 

(D.20) 

(D.21) 


dh 


_ *1 _ 

(1  ♦  h)  sin2(y0) 


(D.22) 


du,/dh  Equation  Solution.  The  Method  of  Separation  of  Variables 
solves  Eq  (D.18). 
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cll 

Solution:  u,  - 

1  1  +  h 


(D.23) 


dat/clh  Equation  Solution.  Rearranging  Eq  (D.19)  gives  the  ODE  a 
form  solved  by  using  an  integration  factor. 


dqt  Qi 


dh  1  +  h 


f1.3Sul.._SLi!L 


{  uo  J  u^  (1  +  h) 


Substituting  known  relationships  for  u0  (Eq  (D.  13)).  q0  (Eq  (D.14)) 


and  u,  (Eq  (D.23))  gives 


dqt  qt  f  3 

dh  *  1  +h  1  2  -  C0JC02(  1  +  h) 


ILL 


C°,2f^il  +  h)  -  C02(l +  h)2 


The  ODE  now  has  the  form  where  an  integrating  factor  is  used  to 
solve  the  ODE.  From  Eq(D.3) 

1  3 


P(h)  - 


Q(h)  -  - 


1 +  h  (1  +  h)  (2  -  C0IC02(1  +  h)) 
C, 


ILL 


Cm2f^~(l  +  h)  -  C02(l  +  h)2 


'01  c 


01 


^2 


From  Eq  (D.3) 


U  •  e  Ip(h,  dh  *  e  jjl  •  h  ■  (1 .  h)  (2 .  Co)Co2(1 .  b)) 


dh 


3 

(2  -  C0IC02(1  +  h))2 
Vl  +  h 
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Thus,  QUOu  -  Q(h)  e 


Jp(h)  dh  gi 


gives 


Q(h)  e 


Jp(h)  dh 


'll 


(2-C01C02(hh))2' 

'/TTh 


C°,2[c^ 1  +  h) "  C02(1  +  h)2j2 

c, 


HL 


"\/c^(l  +  h)2 
Integrating  the  above  equation  gives 

C, 


/qUi)  e  Jp(h) 


dh 


dh  - 


ILL 


Vc^U  +  h) 


Thus,  Eq  (D.4)  gives 


Solution:  q,  - 


«ewVT7h 

yC0l(i  *  h)  12 


(2  -  C01C02(1  +  h))2 


(D.24) 


d^/dh  Equation  Solution.  The  solution  to  this  ODE,  Eq  (D.20),  is 

trivial. 

Solution:  I, -C, 3  (D.25) 

dfy/dh  Equation  Solution.  The  solution  to  this  ODE,  Eq  (D.21),  is 

trivial. 

Solution:  Q,  -  CH  (D.26) 
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daj/dh  Equation  Solution.  This  ODE  is  solved  by  first  rewriting  the 
term  yj  in  Eq  (D.22)  in  terms  of  q(.  From  the  assumed  expansions  for  the 
outer  variables  (Eq  (4.6))  and  using  the  definition  of  q  -  cos(y) 

q  “  q0  +  q l  +  OU2)  -  cos(y)  -  cos(y0  ♦  y,  +  0(e2)) 

The  small  perturbation  expansion  for  the  cosine  function,  Eq  (C.8),  gives 
cos(y0  ♦  y,  +  0(e  ))  -  cos(y0)  -  y,sin(y0)  €  ♦  O(r) 


Equating  the  above  two  expressions  gives 
q0-cos(y0) 


(D.27) 


'  '  sin(y0) 


(D.28) 


Therefore,  Eq  (D.22)  becomes 


dh  (1  ♦  h)  sin3(y0) 


Next,  the  sin3(y0)  term  in  the  above  equation  is  expressed  in  terms  of  q0, 
where,  by  definition,  q0  -  cos(y0). 


s«n3(y0)  2  |  ’  2| 

(1 -cos2(y0))2  (l-q02)2 


But  from  Eq  (D.  14) 

2  1 
Qo  *  ~2 - - 

r~il ♦  h)  -  C02(l  ♦  h)^ 

coi 
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After  combining  the  above  two  equations  and  simplifying,  the 


expression  for  l/sin3(y0)  becomes 


1 


sin3(y0) 


(  2(1  +  h)  -  CftlC02(l  *  h)2 


1”  C0I  -  2(1  -  h)  -Co,Co2(l  -  h)  J 


Substituting  this  relation  into  the  ODE,  Eq  (D.22)  gives 


da, 


di 


dh  (1  +  h) 


2(1  *h)-Cn,Cft2(l  +  h)2  W 


C0I  +  2(1  +  h)  -  C01C02(1  +  hrj 


But  from  Eq  (D.24),  q,  is  known.  Therefore,  after  substitution,  the 
above  ODE  becomes 


1 

(  2(i>h)-c0,cQ2u*hr  v 

(Mi) 

I'Cqi  *  2(I*h)  -  C0tC02(l+h)*J 

^  cn 

VCni(Uh) 


\  "  C01C02(l+h)) 


} 


After  simplifying,  the  ODE  finally  becomes 


-r^*C,2(l  ♦  h) 

_ 

3 

♦  2(l+h)  -  C01C02(l+h)2)2 


The  Method  of  Separation  of  Variables  is  used  to  solve  this  ODE. 
From  Eq  (D.l ),  the  solution  has  the  form 
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CO  I  CM 


al 


_ dh _ 

(-C01  *2(1  +  h)  -  C0jC02(  1  +  h)2) 


(1  +  h)  dh _ 

3 

+  2(1  +  h)  -  C01C02(1  +  h)2)2 


The  first  integral  has  the  following  form  and  solution  (Beyer, 
1984:255),  where  x  -  l  +  h,  dx  -  dh,  a  -  -C0I>  b  -  2  and  c  -  -C01C02. 


2(2cx  +  b) 

(4ac  -  b2Wa  +  bx  +  cx2 


or 


•  _ dh _ 

3 

(-C0I  >2(1  .h)-c0,c02(l  *h)2)* 

( >  -  C0I2C02)  V-C0|  ‘  2(1  *  h)  -  C01C02(1 "  h)2 

The  second  integral  has  the  following  form  and  solution  (Beyer, 
1 984:256 ),.where  x  -  1  +  h,  dx  -  dh,  a  -  -C0|l  b  -  2  and  c  -  -C01C02. 

f  x  dx  2(bx  +  2a) 

- 3  -  - - r  j  2  +  c  or 

/  y\?  (4ac  ’  o  ) va  +  bx  +  cr 

(a  ♦  bi  +  cx  J4 
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(1  ♦  h)dh _ 

3 

(-C0I  .2(1  H-  h)  -  C„,c02(l  ^  h)2)2 

_ (1  *ti)-Cn, 

( I  -  C0,2C02)V-C01  *  2(1  ♦  h)  -  C01C02(1 .  h)2 

Substituting  the  above  two  integrations  into  the  previous  expression 
for  a,  gives 

aj  .  .  _£u _ _ _  -C01C02(1  +  h)  +  1 

'  ^01  ( I  -  C„,2C02 )  V-C0l*2(Wh)-CO|C02(l  .h)2 

_ C|2((l*h)-C0l) _ 

( 1  -  C»|2C02)  V-Cq,  *2(1  ♦  h)  -  C01C02(1  "  h)2  *  15 


Simplifying  the  above  equation  gives 


Solution:  a,  - 


♦C15(D.29) 
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Appendix  E: 

Derivation  of  Solutions  for  Inner  Expansions  of  the  Equations 

of  Motion 


Methods  of  Solution 

The  solutions  to  the  two  sets  of  five,  coupled,  first  order,  linear 
nonhomogeneous  ordinary  differential  equations  (ODEs)  found  in  this 
appendix  are  derived  below.  Five  methods  are  used  to  solve  these 
equations.  The  first  method  used  is  the  Method  of  Separation  of  Variables, 
where  the  ODE  and  its  solution  are  given  as  (Beyer,  1984:315) 

-»  J  f<y>  dy  ♦  J  Hi)  d»  -  K  (H.l) 

K  is  the  constant  of  integration. 

The  second  method  used  to  solve  ODEs  found  in  this  appendix  is  the 
Method  of  Integrating  Factors,  where  the  ODE  has  the  form 

^*P(x)y-Q(i)  (E.2) 

The  integrating  factor  has  the  form 

u  -  e  Jp(,)  41  (E.3) 

and  the  solution  is  given  as  (Beyer,  1984:315) 
ye  dx  -  ]q(x)  e  d*  ♦  K  or 
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The  fourth  method  occurs  when  the  ODE  does  not  have  an 
elementary  function  as  its  solution.  To  derive  a  solution  to  the  ODE, 
judicious  approximations  are  made,  using  binomial  and  Taylor  expansions 
and  assumptions  based  on  the  physics  of  atmospheric  entry  and  on 
experience. 

The  fifth  and  final  method  to  solve  the  following  ODEs  recognizes 
that  some  of  the  ODEs  are  coupled  and  cannot  be  solved  independently. 
Thus,  an  intermediate  solution  is  found  to  two  (or  more)  ODEs  and  by  using 
this  intermediate  solution,  the  original  ODEs  are  solved  (Karasopoulos, 
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1988:177).  For  example,  suppose  there  are  two  ODEs  coupled  and  cannot 
be  solved  independently,  dy/dx  -  f(x,y,z)  and  dz/dx  -  f(x,y,z).  If  the  two 
ODEs  are  combined  by  a  mathematical  operation,  like  addition  or  division, 
the  independent  variable,  x,  may  possibly  disappear,  resulting  in  an  ODE, 
dz/dy  -  f(y,z),  which  is  solvable.  Assuming  there  is  a  solution  to  the  above 
ODE  in  the  form  z  -  Fl(y,z)  or  y  -  F2(y,z),  the  two  original  ODEs  are 
expressed  as  functions  of  only  two  variables  and  thus  are  solvable. 

Constants  of  Integration  Notation 

As  shown  in  Section  IV,  there  are  ten  ODEs  expanded  from  the 
equations  of  motion.  Accompanying  the  solution  of  each  of  these  ODEs  is  a 
constant  of  integration.  To  simplify  the  bookkeeping  involved  in  defining 
all  these  constants,  the  following  notation  is  used  for  the  outer  expansion 
solutions: 

The  constant  of  integration  is  given  as  Kjj,  where  i  denotes  the  order 
(e*  —  €°  or  e1)  of  the  solution  and  j  denotes  the  variable  associated  with  this 
constant. 

j _ variable 

1  u 

2  q 

3  I 

4  Q 

5  a 

Thus  K0)  is  the  constant  of  integration  for  the  variable  u0  and  ICM  is 
the  constant  integration  for  the  variable  In  addition,  the  coupling 
present  in  the  first  order  inner  ODEs  generates  lengthy  constant  terms  in 
their  solutions.  The  following  shorthand  notation  is  used  to  describe  these 
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terms.  Kljk  are  the  constants  multiplying  the  terms  in  the  integral  solution. 
As  before,  K  signifies  an  inner  expansion  solution,  T  signifies  a  first  order 
solution,  j'  signifies  the  variable  being  solved  and  k‘  signifies  the  kth 
constant  being  defined  in  this  solution.  KM,  is  the  first  multiplicative 
constant  for  the  solution  of  variable  Qj.  This  shorthand  will  only  be  used 
for  the  first  order  inner  solutions. 

Inner  Expansion  e°_  Terms 

The  differential  equations  corresponding  to  the  equation  of  motion 
terms  were  derived  in  Section  IV  and  are  repeated  below: 


du0 

-2Bu0e'*(l  +  Atan(y0)) 

(E.8) 

<15  - 

sin(y0) 

d^  “ 

-\Be"< 

(E.9) 

dl0 

B6e~^  cosfao) 

(E.10) 

’ 

sin(y0)  cos(y0) 

d^ 

B6e“*  sin(ao) 

(E.l  1) 

sin(I0)  sin(y0)  cos(y0) 

doo 

BSe"*  sin(ao) 

(E.12) 

^  " 

tan(I0)  sin(y0)  cos(y0) 

Equation  Solution.  The  Method  of  Separation  of  Variables 

solves  Eq  (E.9).  Rearranging  Eq  (E.9)  into  form  of  Eq  (E.l)  gives 

^ -  -  ABe~<  -»  dq0  -  -  \Be_?  d*  -»  J dq0  -  -  J \Be"<  d^ *  K02 
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Solving  the  above  ODE  gives 
Solution:  q0  -  K  Be~^  +  K02 


(E.13) 


duQ/d^  Equation  Solution.  This  ODE  is  solved  by  first  rewriting  the 
terms  sin(y0)  and  cos(y0)  in  Eq  (E.8)  in  terms  of  q0,  where  q0  -  cos(y0). 
Thus.  cos(y0)  *  q0  and  sin(y0)  -  ^Jl-cos2( y0)  -  V 1  -q02.  Substituting  these 
into  the  ODE  gives 


77^  -  -2Bufte~^  [77-7—7  + 


0  (sin(y0)  cos(y0)J 


*  >  -  >-0  1 


-2Bu0e" 


amV  q° 


As  derived  in  Eq  (E.13),  q0  -  q0(^).  To  simplify  the  above  integration, 
the  independent  variable  is  changed  from  l  to  q0.  Since  q0  -  A  Be  *  +  K02* 
dq0«  -i\Be~*  d£  Thus,  the  ODE  now  becomes 


dUfc  2u^  ^  l 


dq0 


M 


Viv  q°, 


Now,  the  Method  of  Separation  of  Variables  is  used  to  solve  the  above 
equation.  Rearranging  the  above  equation  into  form  of  Eq  (E.l )  gives 


^  2u0 1  I  *  ' 

dq„  "  K 


>v  q° 


-I 


du0  2  (  1 


Ur 


yW  q» 


<*q0 


■» 


fdUs 

(2 

(  1  K  \ 

J  Uo’J 

K 

V<V  q°, 

dq0  +  K, 


Oi 


Solving  the  above  integral  (Beyer,  1984:252)  and  remembering,  by 
definition,  (Eq  (D.27))  q0  -  cos(y0)  or  y0  -  Cos' 
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2 

ln(uo)  -  -  fc  Cos  l(q0)  +  2  ln(q0)  +  K0I 


Rearranging  the  above  equation  gives  the  final  solution,  where 
exp( K0l )  -  K01 


Solution:  u0  -  K01  q02  expj-^Cos'^qo) 

=  K0iq02  exp[-|y0] 

-  K01(ABe“?  +  K02)2expj^-^Cos',(ABe"^  K02) 


(E.14) 


qipAK:  Equation  Solution.  The  last  three  ODEs  are  coupled  and  are 
solved  as  such.  First,  the  I0  and  ODEs  are  combined  to  express  I0  as  a 
function  of  c^.  Using  the  chain  rule  for  differentiation 

doo  “  dH;  dcto 


Substituting  for  the  two  known  ODEs  (Eqs  (E.10)  and  (E.12))  on  the 
right  hand  side 

dip  B6e~^  cosfafl)  tan(I0)  sin(yQ)  cos(y0)  tan(IQ) 
dctQ  *  sin(y0)  cos(y0)  B6e‘?  sinto^)  ’  ’  tanfc^) 

Separation  of  variables  solves  the  above  equation.  Rearranging  the 
above  equation  into  form  of  Eq  (E.l )  gives  I0  as  a  function  of  Oq  (Beyer, 
1984:260). 
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dag  tan(op)  tan(I0)  tan(cip) 

.»  [4,. 


-»  ln(sin(I0))  +  ln(sin(a0))  -  ln(sin(K03)) 


■»  sin(I0)  sin(ctp)  -  sin(K03)  or 


,/sin(K02)\  ,/sin(ICfl3) 

■»  - Sin  ( and  «*•  ■  Sin  (iwu 


Next,  the  Qq  and  otp  ODEs  are  combined.  Using  the  chain  rule  for 


differentiation 


dqo  dqo  AL 

dfy  *  dS  d^o 

Substituting  for  the  two  known  ODEs  (Eqs  (E.  11)  and  (E.12))  on  the 


right  hand  side 


dap  B6e  *sin(qp)  sin(Io)sin(y0)  cos(y0) 

dQp  *  '  tan(Ip)  sin(y0)  cos(y0)  B6<f  *  sin(op)  ’  ’  C°S  0 


But  from  Eq  (E.15),  I0  is  known.  Using  this  and  the  trigonometric 
identity  cos(sin'1(x))  -  n  1  -  x2  (Beyer,  1984:141),  the  above  expression  is 


transformed  as  follows. 

-  cos(Ip)  -  -  cosfsi 


/s.-esii 


_  /sin(K03)>z 
l  sinlotf,)  I 


’/a 


sin*(ap)  -  sin2(K03) 


sin^Op) 


'h 


Separation  of  variables  solves  the  above  equation.  Rearranging  the 
above  equation  into  form  of  Eq  (E.l )  gives 


dQp 


sin4(ap)  -  sin2(K03) 


2 

sin  (op) 


►Op 


’/a 


dQp- 


sin(ap)  dap 


Vsin2(otp)  -  sin2(fC03) 


sin(ap)  dap 


Vsin2(op)  -  sin2(K03) 


♦  K, 


04 


Using  the  substitutions  x  -  cos(ap),  dx  -  -  sin(ap)  dap  and  sin2(op) 

-  1  -  cos2(a p)  -  J  -  x2,  the  above  integral  is  easily  solved  (Beyer,  1984:252) 
and  gives  Qp  as  a  function  of  exp. 


Oo- 


‘ \jsin 2 


sin(ap)  d 
(ctp)  -  sin*(K03) 


2r-  *  +  K°4 


dx 


Vl  -i2-  sin2(K03) 


+  K, 


04 


dx 


'  +  Kr 


Vcos2(K03)-  x2  04 


Mss**"" 
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(E.16) 


cos  (K 0-4  -  Op) 


COSlotp) 

cos(K03) 


Finally,  the  dl0/d^  ODE  is  transformed  to  dl0/dy0  via  the  chain  rule 
for  differentiation,  use  of  Eq  (E.9)  and  the  definition  q0  -  cos(y0) 


_  dq„  d(cos<y0)) 

-*Be  -ir~ Tt, — 


dy0 


-  sin<Yo) 


dyQ  ABe  f 
d£  "  sin(y0) 


Using  this  relation  and  Eq  (E.10)  gives 

dl0  dl0  _d£_  B6e~^cos(ap)  sin(y0)  Scoslotp) 
dy0  “  d?;  dy0  "  sin(y0)  cos(y0)  \Be~?  "  Ncos(y0) 


Substituting  the  relation  for  sin(otp)  in  Eq  (E.15)  gives 


dk„_! 


1  - 


sin2(K03)T/2 


sin^(Ip)  J 


dy0  Acos(yp) 


6  [ sin2(I0)  -  sin2dCp3))  >2 
A  sin(I0)  cos(y0) 


Separation  of  variables  solves  the  above  equation.  Rearranging  the 
above  equation  into  form  of  Eq  (E.l)  gives 


^  sec(y0)  dy0  - 


^  J  sec(y0)  dy0  - 


sin(I0)  dl0 


Vsin2(I0)  -  sin2(K03) 
sin(I0)  dip 


or 


Vsin2(I0)  -  sin2dC03) 


05 


The  first  integral  is  easily  solved  (Beyer,  1984:260) 


6 

A 


J  sec(y0)  dy0 
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Using  the  substitutions  x  -  cos(I0),  dx  -  -  sin(I0)  dl0  and  sin2(I0) 

-  1  -  cos2(I0)  -  1  -  x2,  the  second  integral  is  easily  solved  (Beyer,  1984:252) 
and  gives 


£  ln|  tan 


TT  Yo 
\4  2 ) 


sin(I0)  dl0 


Vsin2(I0)  -  sin2(K03) 


-  K, 


05 


dx 


V 1  -  x2-  sin2(K03) 


-  K 


05 


dx 


a/cos2(K03)- 


X 


2  K05 


,/  cos(I0) 

Cos  55s£)J  -  or 


cos(I0)  -  cos(K03)  cos 


f6  .  L  fa  y$\ 

[h  ,nH4 f  2) 


+  K 


05 


Solution:  I0  «  Cos'1  cos(IC03)  cosJ|  In^tanp^  ♦ j  ♦  K05 


Cos 


-1 


cos(IC03)cos|  7  ln|  tan(7  ♦ 


n  Cos’^ABe'^K^ 


+  K 


05 


(E.17) 


dofr/dE  Equation  Solution.  This  solution  is  derived  above  and  is 
repeated  for  continuity.  From  Eqs  (E.15)  and  (E.17)  Oq  is  expressed  in 
terms  of  the  independent  variable,.^. 

,/sinO^h 


Solution:  c*  - 
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repeated  for  continuity.  From  Eqs  (E.16)  and  (E.17)  Qq  is  expressed  in 
terms  of  the  independent  variable..!;. 


Solution:  -  K04  -  Cos"1^^^  -  1CM  (E.19) 


u0ylcos(y0)(l  ♦Man(y0)) 
sin2(y0) 


cos2(y0)  sin(y0). 


(E.20) 

(E.21) 
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dl 

dOj 


-  B6e'* 


yicos(ao)Li,  »-:~x 


aiSinlap) 


BSe 


(cos2(y0)  sin2(y0)J  sin(y0)  cos(y0) 
1  \ 


(E.22) 


yisintap)/  l _ 

[  sin(I0)  [cos2(y0)  sin2(y0)j 


oncost  otp) 


da. 


Ii  sinto^)  cos(I0) 
sin2(I0)  sin(y0)  cos(y0)  sin(I0)  sin(y0)  cos(y0) 

1  .. 


(E.23) 


tan(y0) 


B6e" 


yiSin(oto) 

f  1  1  ) 

tan(I0)  | 

cos2(y0)  sin2(y0)J 

I,  sinCcxp) 


a.cosioo) 


sin2(I0)  sin(y0)  cos(y0)  ^n(I0)  sin(y0)  cos(y0) 


(E.24) 


daj/dE  Equation  Solution.  The  Method  of  Separation  of  Variables 
solves  Eq  (E.21).  Substituting  relations  for  u0  (Eq  (E.13))  and  q0  (Eq  (E.l 4)) 
into  Eq  (E.21 )  and  simplifying  gives 


dq. 

df’qo 


1  ♦ 


SiL 


u 


0) 


/ 


-1  + 


1 

(\Be“< 

*  K02)2 

<0.( 

^Be"^+K02)2expj 

-j-Cos’ 1 1 

W'.K02) 

-y 

(\Be 


-f, 


j^Cos'1  (?\Be~*  +  K02)j 


-  1 


“0! 


Rearranging  the  above  equation  into  form  of  Eq  (E.l )  gives 
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,  f  (xo  -t  „  \[exp  fcos_,(NBe'*  + K02)  ] 

|  dq i  -  J  (AB€  ^K02)\  ^  *  02/J  Jd^K.2 

q,  -  \Be'*  -  K02^ 

♦  ~  |  (\Be-?  +  K02)  expos'1  (\Be~*  ♦  K02)j  d^  *  Kl2 

Using  the  substitutions  x  -  K  Be"^  ♦  K02,  x  -  cos(<j>)  and  the 

.  .  .  .  2sin(<t»)  (to  -  d\  ,(4>  +  d\ 

trigonometric  relationship  cos(d) .  "  COt(  2  ] +  cot[  2  j 

(Karameheti,  1966:624),  the  above  integral  is  reduced  to 
|  (\Be~*  ♦  IC02)  expj^Cos’1  (xBe"?  ♦  IC02)j  d^ 

’  {ic^i  eip['fCos  ,(l)] dx 

d<|>  (cos(d)  -  IC02) 

-If  ,  A  (to  -  d\  fd>  ♦  dYI  (2  \ 

-  ~2  cos(4>)  cotl  —  ♦  cot^ 2~  exP  ^4>j  d4> 

The  above  integral  can  be  reduced  to  a  sum  of  simpler  integrals,  but 
one  resulting  integral  has  an  integrand  of  cot(z)  ebz,  where  z  is  the 
integration  variable,  a  function  of  and  b  is  a  constant.  Using  integration 
by  parts,  the  resulting  integral  has  an  integrand  of  ebz/sin2z,  whose  integral 
cannot  be  expressed  as  a  finite  sum  of  elementary  integrals  (Gradshteyn 
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and  Ryzhik,  1980:197).  Thus,  the  approximation  made  to  solve  the  above 
integrals  is  given  as 

J (  \Be~*  ♦  K02  )exp|j^Cos‘ 1  (  ♦  K02 ) J  d£ 

-  |(\Be"5*  Kojjexp^^Cos'^Koj)  d^ 

Figure  El  compares  the  relationship  between  the  exact  function  (the 
integrand  on  the  left  hand  side  above)  and  the  approximate  function.  While 
there  is  a  significant  disparity  at  smalt  values  of  ^  (low  orbital  altitudes),  the 
above  difference  is  negligible  when  both  functions  are  integrated  from  ^ 

(orbital  altitude  at  the  beginning  of  atmospheric  entry)  to  I;  -  0.  Figure  E2 
shows  the  integrals  of  the  exact  and  approximate  functions  are  identical 
until  5-»0,  where  the  two  functions  deviate  slightly.  This  trend  is  reflected 
in  Figure  E.3,  which  compares  the  exact  solution  of  qlf  obtained  from 
numerical  integration  of  qj,  and  the  approximate  function  derived  for  q  j. 
Thus  substituting  the  approximate  function  into  the  above  solution  for  q, 
gives  an  approximate  solution  for  q,. 

Solution:  q,  -  Km  (k02S  -  hBe"*)  ♦  K,2  (E.25) 

expf^Cos',(K02)| 

where  K,2,  - — ^ — z - ^ -  1  (E.26) 

koi 

dU|/d£  Equation  Solution.  Rearranging  Eq  (E.20)  gives  the  ODE  a 
form  solved  by  using  an  integration  factor. 


135 


Normalized  Functions 


Figure  E 1 .  Comparison  of  Approximate  and  Exact  Functions 
in  Deriving  the  q(  Solution 


du,  (2B (l  +  Man(y0))e  ? 
1  sin(y0) 


u,  -  -u0  ♦ 


2Bu0ylcos(y0)(  1  +  Atan(y0))e' 


sinz(y0) 


2B\  u0  yje~^ 
cos2(y0)  sin(y0) 


The  ODE  now  has  the  form  where  an  integrating  factor  is  used  to 


solve  the  ODE. 


P(tf- 


2B(  1  +  Man(y0))e' 
sin(y0) 
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Normalized  Integrals 


Figure  E2.  Comparison  of  the  Integrations  of  the  Approximate 
and  Exact  Functions  in  Deriving  the  q,  Solution 


Q(^)  -  -u0  * 


2Bu0y1cos(y0)(  1  ♦  Atan(y0))e~* 
sin2(y0) 


2BAu0y1e  * 
cos2(y0)  sin(y0) 


From  Eq  (E.3) 

i>  -  e  Jp(^  -  exp 


f  2B(  1  ♦  Man(y0))e~^ 
sin(y0) 


Factoring  the  trigonometric  terms  and  noting  cos(y0)  -  q0  and  sin(y0) 
V 1  -  qo2.  the  above  integral  in  separated  into  the  sum  of  two 
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Figure  E3.  Comparison  of  the  Approximate  and  Exact  Integral 

Solutions  for  q( 


integrals.  Substituting  in  the  derived  expression  q0  -  K  Be  *  ♦  JC02 
(Eq  (E.I3)),  gives 


u  -  exp 


2Be 


2Bh£l 

1  -  (?\Be~5*K02)2  ABe  ^  ^  I 


Using  the  substitution  x  -  K  Be-*  +  K02,  the  above  integral  is  easily  solved 


v  -  exp 


\  Cos-1  (  \Be-*  ♦  K02  )  -  2  In  (  \Be-*  ♦  K02  ) j 
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(*Be'?«K02)2 


exp 


Thus,  after  factoring  the  trigonometric  terms  and  using  the  definition 
for  q0  as  above,  and  substituting  the  expression  for  qj  (Eq  (E.25)),  where  by 

definition  (Eq  (D.27))  y,  -  -qj/sintyo),  Q(£)u  -  Q(0  e  Ip(^  gives 


exp 


V 


-u0> 


2Bu0y1cos(y0)(  l+\tan(y0))e~*  2B^  u0  y,e'^ 

sin2(y0)  cos2(y0)  sin(y0) 


-K 


2BMC0,e"*  r  1 

+  /  -r  \2  +  K02  KI2l^+  K12 J 

V^Be  '  +  K02) 


oi 


2BK0,e  ^(^*Be  ^  +  K02)  r  ,  1 

3“  [-K12|\Be  +  K02  K12,^  +  Kl2] 


1  -  (\Be"^K02)2p 


Distributing  the  terms  in  the  above  equation  and  integrating  gives 
(where  Aj  -  A6  are  constants) 

(\B?-<.K0j)2 


|  Q(0«  d?  -  K„  -  K„,5  ♦  A, 
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The  first  and  second  integrals  above  are  easily  solved  by  using  the 
substitution  x  «  \  Be"*  +  K02.  The  third  integral  is  solved  by  using  the 

following  relation  (integration  by  parts  and  (Beyer,  1984:278)). 


The  fourth  and  fifth  integrals  are  easily  solved  using  the  substitution 
x  -  A  Be"*  ♦  K02  and  the  following  relation  (Beyer,  1984:253,254). 
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The  sixth  integral  is  solved  by  using  the  substitution 
x  -  and  the  following  relation  (integration  by  parts  and  (Beyer, 
1984:257)). 


I 

I 

I 

I 


Thus,  Eq  (E.4)  gives  the  solution  as  Uj  - 


fo(i)  e  H  « >  K 


fP(0  K 


or 
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Solution:  u(  -  (A.Be_i  *  K02 )2  expf-  ^ Cos'1  ( <\Be"' »  K02) 


K  ,g^0l(^l2^02^ljl)iil  11  2K|2|  K0|  ,/  ,  \ 

JtBe-^K.,  "  *  Sm  ,  Be  K“" 


f|K0,K,2-K0,K12|](K025^Be-<)  2»  ,,  K 


1  -  (\Be"*  +  K02)2 


AVI  -Kr 


2 V ( 1  -K022)  [  1  ■ -ft Be"<  ♦  K02)2  ]  -  2K02  (  ABe~*  +  K02  )  ♦  2  . 

x  In  - —7 -  (E.27) 

L  ABe  *  J 


ctaj/dE  Equation  Solution.  The  last  three  ODEs  are  coupled  and  are 
solved  as  such.  First,  the  a,  ODE  is  solved  using  judicious  assumptions. 

This  solution  is  used  to  derive  an  expression  for  Ij.  Finally  the  solutions  for 
a,  and  Ij  are  used  to  derive  an  expression  for  fy.  The  expression  for  da,/d^ 
is  given  as  Eq  (E.24). 


“  cot(y0)  - B6c  tan(I0)  (cos2(yo)  -  sin2(yo) 


I,  sin(Oft) 


a,cos(an) 


sin2(I0)  sin(y0)  cos(y0)  tan(I0)  sin(y0)  cos(y0) 


Since  re-entry  vehicles  enter  a  planetary  atmosphere  with  small  bank 
angles  (a),  the  value  for  6,  where  by  definition  (Eq  (3.33))  6  -  (Ci/Cp) 

sin(o),  is  very  small,  or  6  -  0(0.001).  Combining  this  with  the  fact  e’^  is  also 
small,  or  e‘*  -  0(0.1),  results  in  a  simplifying  assumption  for  the  da,/d(; 

ODE.  Assuming  cot(y0)  »  B6e*l  which  is  valid  for  y0  <  80*  and  typical 


values  for  B  (B  -  0(10))  (Hillje,  1969:2),  the  coupling  in  the  da,/d(;  ODE  is 
removed.  Thus,  the  da,/d^  ODE  becomes 


cot(y0) 


cos(y0) 

sin(y0) 


This  ODE  is  solved  by  first  rewriting  the  terms  sin(y0)  and  cos(y0)  in 
Eq  (E.8)  in  terms  of  j,  where,  q  -  cos(y0)  -  \  Be'*  ♦  IC02.  Thus,  sin(y0)  - 
Vl-cos2(y0)  -  A/HiNBe'*  +  K02)2.  Substituting  these  into  the  above  ODE 
gives 


da,  \Be~*  »  K02 
d*  *  K02)2 


Now,  the  Method  of  Separation  of  Variables  is  used  to  solve  the  above 
equation.  Rearranging  the  above  equation  into  form  of  Eq  (E.l)  gives 


a 


I 


[  <\Be~*  ♦  K02 

I 

j  V  l-^Be'**K02)2 


Using  the  substitution  x  -  ^Be-^,  the  above  integral  is  reduced  to 


f  dX  K 

f  dX 

VU  -  K022)  -  2K02I  - 12 

W(  1  -  K022)  -  2K02x  -  x2 

The  first  integral  is  easily  solved  (Beyer,  1984:255)  and  gives 


r  di 

.  V<1  -  K022)  -  2  K02i  -  xz 


With  some  manipulation,  the  second  integral  is  easily  solved  (Beyer, 
1984:257)  and  gives 
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Substituting  the  above  two  integrals  back  into  the  original  expression 
for  a,  gives 


Figure  E4  compares  the  exact  solution  of  aj,  obtained  from  numerical 
integration  of  alt  and  the  approximate  function  derived  for  aj.  While  there 
is  a  minor  disparity  at  small  values  of  I;  (low  orbital  altitudes),  the  above 
difference  is  negligible  when  the  magnitude  of  the  disparity  is  compared  to 
the  magnitude  of  the  functions  at  this  altitude.  This  difference  between  the 
exact  and  approximate  solutions  reflect  the  same  behavior  as  found  in  the 
approximate  solution  for  q,  (Figure  E.2) 


Normalized  Solutions 


Figure  E4.  Comparison  of  the  Approximate  and  Exact  Integral 

Solutions  for  a. 


The  above  ODE  is  transformed  to  dl,/dy0  via  the  chain  rule  for 
differentiation,  use  of  Eq  (E.9)  and  q0  -  cos(y0) 


do.  d(cos(y0)) 


_  _ f  «mo  -v— ov  .  ,  ,  dy0 

-*Be  ■'df • — 3? — *-s'n(yo>-df  - 


dy0  \Be  ^ 
d^  "  sin(y0) 


Using  this  relation  and  Eq  (E.22)  gives 

dl |  dl,  d^ 
dy0  *  d^  dy0 


sin(y0) 


B6e”* 


yjeosCao) 


/  1  1  > 

a^intoto) 

kcos2(y0)  sin2(y0); 

’  sin(y0)  cos(y0) 
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&[  ,  .  sin(yQ)  l  ^  aisiniap)' 

M\  Yicoso^)  cos2( y  ) "  sin(y 0)  "  cos(y0) 

From  Eqs  (E.18),  (E.25)  and  (E.28),  relations  for  Oq,  qt  and  a,  are 
given,  and  by  definition  (Eq  (D.28)),  y,  -  -qj/sin(y0).  Substituting  these 
relations  into  the  above  ODE  gives 


Before  separation  of  variables  is  used  to  solve  the  above  equation, 
several  approximations  are  used  to  bring  the  integration  task  down  to  a 
manageable  level.  First,  using  order  of  magnitude  analysis  involving  6/\, 
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and  K05,  6A  -  0(0.00 1 )  and  K05  >  6A .  Therefore,  the  following 
approximation  is  made. 

sin  |  ln|tan|j +  )j +  *05  ^  sin(K05) 

Therefore 

cosj^  ln|tan^  + +  *05  ~cos(K05) 

Also,  through  numerical  analysis,  the  following  approximations  are  made. 
ln(cos(y0)  -  K02)  -  -  cos(y0) 


Substituting  these  relations  into  the  above  ODE  and  combining 
integrands  gives 

k{[  kI2I  (*020  *  In(^B))  -  (1  ♦  K02)cos(y0))  ♦  KI2j 


X  /KQ2  2  cos(yo)  +  KI5  *  Yo  •  f 

.yl-Ko2  1 


sinOC03) 

■  i  ...rr.i.rim.  ■ ,  ",  afeBsaggas  ..ju-ji 

cos(y0)  V 1  •  sin2(K05)  cos2(K03) 
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Rearranging  the  above  equation  into  form  of  Eq  (E.l )  gives 


h-*i3  +  A 


Kq2  sin(KQ3) 


*^02  j* 

V  i-K022  V 1  -  sin2(K05)  cos2(K03)  Y° 


sin(lC05)  1C 1 21  (Kq2  +  l )  cos(K03) 
V 1  -sin2(K05)cos2(K03) 


f  cos(y0) 
sin2(y0)  Yo 


sin(K05) 

Ki2+K02Ki2i  ( l+ln(i\B)) 

]  cos(K03) 

r  dYo 

V 1  -  sin2(K05)  cos2(K03) 

j 

sin2(y0) 

sin(IC05) 

^12+^02  ^121  (l+ltl(^B)) 

]  cos(K03) 

r  d^o 

Vl  -  sin2(K05)  cos2(K03) 

1  j 

cos2(y0) 

sin(K0S)  1C|2|  cosOCqj)  (Kq2  +  l)  +  sin(K03) 


(JL  v  ) 

2  '  N5 


Vl  -  sin2(K05)  cos2(K03) 


dy0 

cos(y0) 


sin(K03) 


_ ==  Yo  dy0 

Vl  -  sin2(K05)  cos2(K03)  J  cos(yo)  j 


The  first  five  integrals  are  easily  solved  (Beyer.  1984:260,261).  The 
sixth  integral  is  given  by  (Gradshteyn  and  Ryzhik,  1980:189) 

oo 

[jisu-.y  iE2t[x2ii,°-1 
J  cos(x)  ^  (2k+n  +1)  (2k!) 

k-0 
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The  term  |E2k|  is  the  Euler  constant:  its  first  few  terms  are  given  as 
(Gradshteyn  and  Ryzhik,  1980:xxix)  |E0|  -  1,  IEjI  -  1  and  |E4|  -  5.  Thus  the 

above  integral  becomes 


Ypdy,  y,*  5yl  8 
cos(y0)  2  8  144  °ly°  ’  2  8 


Substituting  y0  -  Cos‘1(\Be  *  ♦  K02)  gives  the  solution  for  I,  as 


i'  fir  *133  ) 

- T - 

+ 

l 

*131  [*132-Cos(y0)j 

sin(y0) 

♦  K 


134 


*02  Yo  Yo2  Yq41 
2  2  8 


“02 


*13IK132 


„  fl r 

+  ”‘134  2  '  R 


IS 


id  tan 


tt  yo) 
4  *  2 


Solution:  I,  -  K13  +  ^ 


KI3I 

L  *133 

l32'\Be-?.K02 

^1  -  ^Be_5  +  K02)2 

- 

n  Cos''1 

/Be-<.K02)f 

*131*132  *  ^134 

. 

[*  '*) 

In 

tan 

k 

2  Jj 

♦  K 


134 


K0J  Cos'1  (  M3e'* .  K02 )  [cos’1  (  \Be‘?  •  K02  )]J 


VT 


L02 
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Cos'^KBe"5*  K02)  4 

■ 

OO 

sin(K05)  cos(KQ3) 

1  -  sin2(K05)  cos2(K03) 


ki32  "  Kt2i  +  ^02) 

^(33  *  ^12  +  ^02  ^121  (  1  +lfl(?sB) ) 

y  _ sin(K03.j _ 

kl34  I  y  "  y 

V 1  -  sinz(K05)  cosz(K03) 


(E.29) 


(E.30) 


Figure  E5  compares  the  relationship  between  the  exact  solution  of  I,, 
derived  by  numerical  integration,  and  the  approximate,  analytical  solution 
derived  above.  As  in  Figure  E3,  the  correlation  between  exact  and 
approximate  solutions  is  good,  with  a  small  difference  appearing  as  £~>0. 
This  is  the  same  phenomena  found  in  the  approximation  for  q,.  Thus,  even 

with  the  three  approximations  used  to  derive  this  solution,  a  significant 
portion  of  the  disparity  between  the  exact  and  approximate  solutions  are 
attributable  to  the  approximate  expression  used  for  q,. 


dfl|/d£  Equation  Solution.  The  expression  for  dQ,/d^  is  given  as  Eq 

(E.23). 


dQ,  , 

’yjsinlao) 

f  1  1  ) 

d^  "  B6e 

sin(I0)  1 

^cos2(y0)  sin2(y0)/ 

li  sinlop)  cos(Iq)  ajCOsCop) 

sin2(I0)  sin(y0)  cos(y0) +  sin(I0)  sin(y0)  cos(y0) 
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Normalized  Solutions 


Figure  E5.  Comparison  of  the  Approximate  and  Exact  Integral 

Solutions  for  I, 


The  above  ODE  is  transformed  to  dfi,/dy0  via  the  chain  rule  for 
differentiation,  use  of  Eq  (E.9)  and  q0  -  cos(y0) 


dq0  d(cos<y0)) 


_ _ f  uwo  jov  .  ,  x  dy0  dy0  \Be~* 

-NBe  -df- — 5? — --sin<yo)-df  -» 


Using  this  relation  and  Eq  (E.23)  gives 


dQt  dQt  (j^ 

dy0 "  tt,  <Jy0 
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+ 


1 

sin(l0)  cos(y0) 


|ajCOs(ag) 


I,  sin(aft)  cos(I0) 
sin(l0) 


From  Eqs  (E.18).  (E.25),  (E.28)  and  (E.29),  relations  for  Oq,  q,.  a,  and 
I|  are  given,  and  by  definition  (Eq  (D.28)),  yj  -  -qj/sin(y0).  Substituting 
these  relations  into  the  above  ODE  gives 


dQ,  dQ, 
dy0“  dJ;  dy0 


6 

y,  sin(y0)sin(cto) 

t  1  i  \ 

X 

sin(I0) 

(cos2(y0)  sin2(y0)J 

1 

+  sin(I0)  cos(y0) 


^coslc^) 


I,  sinlap)  cosOp) 
sin(I0) 


K|33"^02^l2l*n(COS^0^"^02)'^l2l  cos^yo^] 


fsin(K03)  sin(K03) 


sin*(yft)  cos*(y0)j 


6 


o 

<M 

C 

vi 

1  - 

1  - 

\ 

cos(K03)  cos 

iln 

) 

|+^05 

■2 

/ 

+ 


sin(K03)  cos(K03)  cos 

i,n 

tan 

\ 

n  Cos‘' 
4  + 

(\Bc-?.K02) 

2 

\\ 

’  ^05 

±1 _ _± 

cos(y0) 

1  - 

cos(K03 

|)  cos 

(|ln(tan( 

4* 

T 

3 

2 

X 


KI33  ) 

Kl3  +  \ 

k»3 1  (^132  -  COS(y0)J 

.  V 

KQ2  YO 

2  4" 

yjl  Yo 

sin(y0) 

V 

+  kl34 

1 - 

1 

© 

ls> 

2  ■  8 

id  tan 


fe.&l 


V 


4  2 


Before  separation  of  variables  is  used  to  solve  the  above  equation, 
several  approximations  are  used  to  bring  the  integration  task  down  to  a 
manageable  level.  First,  using  numerical  analysis  involving  6A  and  Kos, 
6/\  -  0(0.001 )  and  K05  »  6M .  Therefore,  the  following  approximation  is 
made. 


sin(K05) 


Therefore 


cos(K05) 
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Also,  through  numerical  analysis,  the  following  approximations  are  made. 
ln(cos(y0)  -  K02)  -  -  cos(y0) 

-  In  (2^/1  -  K022  sin(y0)  -2  K02  cos(y0)  +  2  )  -  0 

Vi-k02z 

Substituting  these  relations  into  the  above  ODE  .combining 
integrands  and  rearranging  the  result  into  the  form  of  Eq  (E.l)  gives  (where 
A |  -  An  are  constants) 


r  d^o  +A 

(■  dy0  r  cos(y0)  dy0 

Al 

j 

sin2(y0)  * 

cos2(y0)  3J  sin(y0) 

f  <*Vg  . 

r  dv0  . 

r  dY0 

<J  cos(y0)  Asj 

cos(y0)  sin(y0)  6 

sin(y0)  cos2(y0) 

7J  cos(y0) 

fyfiidyg  f  yQ4JyA  r 

cos(y0)  A*J  cos(y0)  AioJ  d*o 

dy0 


The  first  six  integrals  are  easily  solved  (Beyer,  1984:260,261,263). 
The  seventh  through  the  ninth  integrals  are  given  by  (Gradshteyn  and 
Ryzhik,  1980:189) 

00 

di  y  jEajif^ 

cos(x)  ^  (2k*n  *1)  (2k!) 

k«0 


+  A 


It 


Id 


11 


cos(y0) 
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The  term  |E21c|  is  the  Euler  constant:  its  first  few  terms  are  given  as 
(Gradshteyn  and  Ryzhik,  1980:xxix)  |E0|  -  1,  |e2|  -  1  and  |e4|  -  5.  Thus  the 

above  integrals  becomes 


ry0n  dYo  _  ypn+l 
cos(y0)  “  n  +  1 


The  eleventh  integral  is  transformed  to  put  the  equation  in  a  suitable 
form.  By  using  trigonometric  identities  (Beyer,  1984:138) 


tan 


TT  y$\  cos(y0) 

4  +  2  J  "  1  -  sin(y0) 


Using  the  law  of  logarithms 


Iru  tan 


( tt  y o'! 

U  *  2  JJ 


cos(y0) 

i.n 

‘1  ♦  sin(y0)' 

.1  -  sin(y0) 

"  2  n 

.1  -  sin(y0)_ 

Using  the  series  expansion  for  the  right  hand  side  term  yields  (Beyer, 
1984:297) 


l“ln| 
2 


1  ♦  sin(yp) 


[1  -  sin(yo)J 


-  £  sin2°''(y»)  -  sin(yo)  *  .  o(sin5(yo)) 


n»l 


2n  -1 


Therefore 


Ira  tan 


ki!  2h 


dy0 


cos(y0) 


r  1  sin3(y0)' 

tan(y0)  *  3cos(y0)J  d*> 


4ln(cos(y0))  cos2(y0) 
3  ‘  6 
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Solving  the  above  integrals  and  substituting  y0  -  Cos'^ABe'*  ♦  K02) 


gives  the  solution  for  I,  as 


v  v  -\ii 
cJ  k141  k132  ~  COs(va)  M  Y. 

♦  -  _ L  _COS{yo’}  *• 

1  14  X  l  sin(vJ  cos(v 


”  _  —  l 

132  -COS(yoV  +  KH2  K133  I  r  fn  Y& 
sin(y0)  cos(y0)  K^3  ln[tan[^4  +  2 


+  (*.32  *H2  +  *144  (f  -  *|5))  [lMgpS(y0))  _  COS^ 


1  -Kr 


V  V  133 

^141  RI32'XB  -f  v 

Solution:  Q,  -KH+f *  — — .  °2/ 4  -  1.12  1 

'  M  A  /  /  _  \2  XRo~s  ^ 


1  -  f *Bc  *+K< 


+  *m  ln[tanl^  + 


n  Cos'^ABe'5*  K02) 


+  *132  *^142  +  *144  2  "  *15 


4  ln(ABe"*  +  K02)  (ABe"*  +  K02)2 


*02  ^144  ^0S  1  (  ^  +  *02  ) 


l  "  Kr 
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+  *"145 

Cos'1 

^Be"^K02)]3  3 

Cos'^Be'S  +  KoJ]5 

6 

40 

*  *"H6 

l - 1 

0 

0 

c/> 

1 

'  NBe"*  +  K02  )]2  [cos'1  (  ?\Be~*  ♦  K02  )]4 

2 

8  JJ 

where  K 


m 


141  sin(K03) 

6  sin(K05)  cos(K05)  sin(Kfta)  cos2(Kft,) 


142 


L03' 


L03j 


Q21 


h  ( 1  -  sin2(K05)  cos2(K03))2 

*M43  •  *"142  (^133  ■  ^132) 

K|3^cos(K03)(sindCQ5)  -  K022) 

1  -  sin2(K05)  cos2(K03) 


sindCoj)  cos(K03)  Kj5  -  2I 
1  -  sin2(K05)  cos2(Kn,) 


03' 


sin(K05)  cos(K03) 

144  1  -  sin2(Kn«)  cos2(Kn,) 


"05' 

K142  tan(K03) 


kl42 

*M5'  sin(K05) 


‘146  r"144 


*"02K143 


VT 


"02 


(E.31) 


(E.32) 
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Figure  E6.  Comparison  of  the  Approximate  and  Exact  Integral 

Solutions  for  Q, 

Figure  E6  compares  the  relationship  between  the  exact  solution  of  Qt, 
derived  by  numerical  integration,  and  the  approximate,  analytical  solution 
derived  above.  As  in  Figure  E3,  the  correlation  between  exact  and 
approximate  solutions  is  good,  with  a  small  difference  appearing  as  £;-»0. 
This  is  the  same  phenomena  found  in  the  approximation  for  q,.  Thus,  even 

with  the  three  approximations  used  to  derive  this  solution,  a  significant 
portion  of  the  disparity  between  the  exact  and  approximate  solutions  are 
attributable  to  the  approximate  expression  used  for  q,. 
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Appendix  F: 

Sample  Earth  Atmospheric  Entry  Trajectories 


This  appendix  presents  an  example  Earth  atmospheric  trajectory  to 
demonstrate  the  use  of  the  zero  and  first  order  composite  expansions 
derived  in  Section  IV,  Eqs  (4.107M4.1 1 1).  The  initial  conditions  are  for  an 
Apollo-type  reentry  vehicle  and  are  derived  from  theoretical  and  flight  test 
data  (Hillje,  1969:2-10) 

For  the  following  sample  trajectory  presented  in  Figures  FI -F5,  the 
independent  and  dependent  variables  are  defined  in  Sections  II  and  III  and 
are  repeated  below. 

y  -  Vehicle  Altitude  (m) 

u  -  Non-dimensional  Speed  Ratio,  u  s  V2  cos2(y)/gsrs 
y  -  Flight  Path  Angle  (deg) 

I  -  Inclination  Angle  (deg) 

Q  -  Longitude  of  the  Ascending  Node  (deg) 
a  »  Latitude  at  Epoch  (deg) 
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Figure  F 1 .  Comparison  or  the  Numerical  and  Analytical 
Solutions  for  the  Speed  Ratio,  u 


I 
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j 


I 

I 

I 


Figure  F2.  Comparison  of  the  Numerical  and  Analytical 
Solutions  for  the  Modified  Flight  Path  Angle,  q 
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I 


Figure  F3.  Comparison  of  the  Numerical  and  Analytical 
Solutions  for  the  Inclination  Angle,  I 


Figure  F4.  Comparison  of  the  Numerical  and  Analytical 
Solutions  for  the  Longitude  of  the  Ascending  Node,  Q 
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LATITUDE  AT  EPOCH 


ALTITUDE  (km) 


Figure  F5.  Comparison  of  the  Numerical  and  Analytical 
Solutions  for  the  Latitude  of  Epoch,  a 
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